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A. Analysis of Agronomic and Spectral Data For Physical Understanding 

Crop identification and area estimation promises to be one of the 
major applications of remote sensing and the Large Area Crop Inventory 
Experiment (LACIE) has pushed the technology to near operational use 
for wheat (MacDonald and Hall, 1978). Remote sensing also offers great 
potential for obtaining accurate and timely information about the 
condition and yield of crops (Bauer, 1975) . 

To fully realize the potential of remote sensing for crop identi- 
fication, condition assessment, and yield prediction it is important to 
understand and quantify the relation of agronomic characteristics of 
crops to their raultispectral reflectance properties. For example, it 
is essential to know in which regions of the spectrum information relating 
to variations in crop parameters is contained. This information is 
necessary for the optimum use of current Landsat technology, as well 
as for the design and development of future remote sensing systems. 

Differences among crop species and dynamic changes due to growth, 
development, stress, and varying cultural practices cause differences in 
the reflectance spectra of crops. Many of the factors affecting the 
reflectance properties of plant leaves have been identified and investi- 
gated utilizing laboratory measurements. The relationships of physical- 
biological parameters such as chlorophyll concentration, water content and 
leaf morphology to reflectance, transmittance and absorption have been 
well-established for leaves. Some of the papers and reviews describing these 
relationships include: Gates et al. (1965), Breece and Holmes (1971), 

Gausman et al. (1970), and Sinclair et al. (1971). 


This section describing the results of Task lA, Analysis of Spectral Data 
for Physical Understanding, was prepared by the task leaders. Dr. Marvin Bauer 
and Marilyn Hixson. Drs. Craig Daughtry, Virgil Anderson, and Vern Vanderbilt 
made valuable contributions to the experiment design, data analysis, and 
interpretation of results. Graduate assistants, John Ahlrichs, Don 
Crecelius., Jeff Kollenkark, and Christina Stellon conducted many of the 
analyses . 


Knowledge of the reflectance characterise ics of single leaves is 
basic to understanding the reflectance properties of crop canopies in 
the field, but cannot be applied directly since there are significant 
differences in the spectra of single leaves and canopies. The reflec- 
tance characteristics of canopies are considerably more complex than 
those of single leaves because there are many more interacting vari- 
ables in canopies. Some of the more important agronomic parameters 
influencing the reflectance of field-grown canopies are: leaf area 

index, biomass, leaf angle, soil cover percentage, soil color, and 
leaf color. Differences in these parameters are caused by variations 
in many cultural and environmental factors, including planting date, 
cultivar, seeding rate, fertilization, soil moisture, and temperature. 
Solar elevation and azimuth angle and the view angle and direction of 
the sensor also affect the measured reflectance of crops and soils. 

The spectral and agronomic measurements which have been acquired 
during the three years of the LACIE field research program are being 
analyzed to provide an understanding of the relationship of reflectance 
to the biological and physical characteristics of crops and soils. The 
primary data being analyzed are the spectrometer data acquired by the 
truck- and helicopter-borne systems. These data are particularly useful 
because the spectral data were acquired in 0.01 ym wavelength intervals 
and are calibrated in terms of bidirectional reflectance faction. Having 
the entire spectrum from 0,4 to 2.4 ym permits simulation of the response 
in any specified waveband. In other words, the analysis is not restrict- 
ed to a fixed set of bands such as Landsat MSS or one of the aircraft 
scanner systems. Calibration of the data permits valid comparisons to 
be made among different dates, locations , and sensors. 


The overall objective of the analyses being conducted is to quanti- 
tatively determine the spectral-temporal characteristics of wheat, 
small grains, and other agricultural crops. The specific objectives 
of the analyses reported here are: 

1. To determine the relationship of agronomic variables 
such as biomass and leaf area index to the multispectral 
reflectance of spring wheat and evaluate the potential 
for predicting the agronomic characteristics of wheat 
canopies from reflectance measurements. 

2. To examine the effects of cultural and environmental 
factors on the spectral response of spring wheat. 

3. To determine the dlscriminability of wheat and other 
crops as a function of acquisition date and spectral 
band selection. 
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PREDICTION OF AGRONOMIC CHARACTERISTICS OF SPRING RTiEAT 
CANOPIES FROM REFLECTANCE MEASUREMENTS 


One of the major long term goals of agricultural remote sensing 
research is to estimate or predict crop variables that can subsequently 
be used to assess crop growth and vigor or be entered into a yield 
production model. To achieve this goal the relationship between the 
spectral reflectance of crop canopies and their biological and physical 
characteristics must be understood. Therefore, the relationships of 
canopy variables to multispectral reflectance were first Investigated. 
This information was then used to develop regression equations for 
predicting canopy variables from reflectance data.. The analyses 
were performed using wavelength bands of the current and future 
satellite multispectral scanner systems. 

1.1 Experimental Approach 

The data used for these analyses were collected using a truck- 

mounted spectrometer during the first two years of the LACIE field 

research program at the agriculture experiment station at Willis ton. 

North Dakota. Measurements were made at approximately weekly intervals 

of the bidirectional reflectance factor and agronomic characteristics 

of 32 plots of spring wheat (Bauer et al. , 1978). The design of 

4 

the experiment was a 2 factorial wj th two replications. The factors 
and levels were : 

Soil Moisture:. (1) wheat during previous year, (2) fallow during 

previous year. 

Planting Date: (1) early, (2) late 

Cultivar: (1) standard height, (2) semi-dwarf. 

Nitrogen Fertilization: (1) none, (2) 34kg/ha 



The factors and levels were selected to represent regional agri- 
cultural practices that affect the growth, development, and yield of 
spring wheat. The treatments resulted in a relatively wide range of 
types of wheat canopies, differing in maturity, biomass, and percent 
soil cover at any one time and over the season. 

In addition to reflectance, the following agronomic variables 
were measured: maturity stage, height, percent soil cover, leaf area 

index, percent green leaves, fresh and dry biomass, and plant water 
content (difference between fresh and dry biomass). Vertical and 
oblique photographs were taken of each plot on each measurement date. 

After processing, graphs of the spectral data were examined to 
verify data quality and qualitatively assess the information contained 
in the spectra. Correlation and regression analyses were used to 
relate biological and physical variables describing the canopies to 
spectral response. 

Since the application of remotely sensed spectral measurements 
will be with multispectral scanner systems which measure the spectral 
response in selected wavelength bands, the statistical analyses were 
performed using wavelength bands. Two different sets of bands were 
considered, Landsat MSS and Landsat thematic mapper. The thematic 
mapper is a second generation satellite multispectral scanner system 
planned for Landsat-D which is to be launched in 1981. The Landsat 
MSS bands are: 0.5-0. 6 (green), 0.6-0. 7 (red), 0,7^G.8(near infrared), 
and 0. 8-1.1 pm (near infrared). The thematic mapper bands which are 
narrower and sample more parts of the spectrum are: 0.45-0.52(blue) , 

0.52-0.60 (green), 0.63-0. 69(red), 0.76-0.90 (near infrared) , 1.55-1.75 
(middle infrared), and 2.08-2. 35pm (middle Infrared). The thematic 
mapper will also have a thermal infrared band, but thermal measurements 
were not acquired for this experiment. 
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1.2 Results and Discussion 

Relation of Canopy to Variables to Multlspectral Reflectance 

The amount of vegetation present is one of the principal factors 
influencing the reflectance of crop canopies. Figure A-1 illustrates 
the effect of amount of vegetation as measured by leaf area index, 
biomass, and percent soil cover on the spectral response during the 
period between tillering and the beginning of heading when the maximum 
green leaf area was reached. As leaf area and biomass increase there 
is a progressive and characteristic decrease in the reflectance of the 
chlorophyll absorption region, increase in the near infrared reflectance, 
and decrease in the middle infrared reflectance. 

The relationships between several agronomicaly important 
characteristics of crop canopies and reflectances in selected wavelength 
bands are illustrated in Figures A-2 to A-5 . These plots include 
data from all treatments when green leaves were present (seedling 
through flowering stages of maturity). Plotting the spectral reflectance 
versus canopy variables allows visual assessment of the relation of 
canopy variables and reflectance, and permits better interpretation of 
the correlation analyses (Table A-1). It should be noted that a 
portion of the scatter in the data is associated with various agronomic 
treatments, as well as measurement errors in the independent and dependent 
variables . 

The relationships between the 0.63-0.69 and 0. 6-0.7 pm reflectances 
and green vegetation results from the strong absorption of incident 
radiation by the chlorophylls . Reflectance in these wavelength bands 
is inversely related to the amount of chlorophyll present in the 
vegetative canopy and thus is sensitive to green or photosynthetically 
active vegetation. 

The relationship of the near infrared wavelength bands , 0.76-0.90 
and 0.8-1. 1pm, results from the lack of absorption and high degree of 
scattering of near infrared radiation by green leaves. Over the ranges 
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Figure A-1. Effect i ofi leaf area index, percent soil cover, dry biomass, 
and plant height on the spectral reflectance of spring wheat 
during the period between tillering and the beginning of 
heading, when the maxinmm green leaf area is reached. Data 
were acquired at Willis ton, North Dakota on May 28 - June 18, 
1976 and Include plots with different soil moisture levels, 
planting dates, nitrogen fertilization and cultivars. 
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Figure A-2. Relationship of leaf area index to reflectances in the 
chlorophyll absorption (0.63-0.69 pm) and near infrared 
(0. 76-0. 90 pm) regions. Measurements for seedling through 
flowering stages of maturity are included for 16 treatments 
representing different levels of soil moisture availability, 
planting dates, nitrogen fertilization, and cultivar. 


Bidirectional Reflectance Factor (%) 



£ 


50r 0.78 -0.90|Jm 


40. i ^ * 



Percent Soil Cover 


Figure A-3. Relationship of percent soil cover to near infrared 
reflectance of spring wheat canopies. 
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Fresh Biomass (g/m^) 

Figure A-4. Relationship of fresh biomass to middle infrared 
reflectance of spring wheat canopies. 
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% Plant Moisture Plant Water Content (g/m^) 


Figure A- 5. Relationship of two measures of plant moisture, percent plant 

moisture and plant water content, to middle infrared reflectance 
of spring wheat canopies. 
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of our data there is a nearly linear relationship between the near 
infrared reflectance and the amount of green vegetation present . 

Reflectances in the middle infrared wavelength bands, 1.55-1.75 
and 2.08-2.35 pm, are strongly influenced by plant moisture. The 
reflectances are inversely related in a nonlinear manner to the amount 
of plant moisture present in the canopy. 

AS a way of summarizing these relationships the linear correlations 
of five different canopy variables x<rith the Landsat MSS and thematic 
mapper wavelength bands are shown in Table A-1. A near infrared 
wavelength band, 0.76-0.90 pm, is most highly correlated with percent 
soil cover and leaf area index. I^ile the middle infrared band, 

2.08-2.35 pm, has the highest correlation of the six individual 
thematic mapper bands with fresh biomass, dry biomass, aiid plant 
water, and is highly correlated with leaf area index and percent soil 
cover. These results Illustrate the importance of making measurements 
in several regions of the spectrum, especially in the middle infrared 
wavelength region which is not measured by the present Landsat multispectral 
scanners. 

The Correlations of the Landsat MSS and the proposed thematic 
mapper bands with the crop canopy variables can be compared in Table 
A-1. The correlation of canopy variables with each Landsat MSS band 
is less than with the corresponding thematic mapper band. Although 
the differences in correlation coefficients between the Landsat MSS and 
thematic mapper bands are relatively small, the coefficients for the 
thematic mapper bands are consistently higher. The higher correlations 
for the thematic mapper bands are attributed to more optimal width and 
location of bands with respect to the spectral characteristics of 
vegetation. 

For example, the data in Table A-1 demonstrate a disadvantage of 
collecting data in the 0. 7-0. 8 pm wavelength region. The inclusion of 
the region of rapid transition from the chlorophyll absorption region 
of the spectrum to the highly reflecting near infrared region (0.70-0.74 pm) 
results in a weaker relation between reflectance and the crop canopy 
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Table A-1. The linear correlations of reflectances in the thematic 
mapper and Landsat MSS wavelength bands with percent 
soil cover, leaf area index, fresh and dry biomass, and 
plant water content. 


Wavelength 
Band (pm) 

Percent 

Soil 

Cover 

1 

i 

Leaf 

Area 

Index 

Fresh 

Biomass 

Dry 

Biomass 

Plant 

Water 

Content 

Thematic Mapper 

0.45-0.52 

-0.82 

-0.79 

-0.75 

-0.69 

-0.76 

0.52-0.60 

-0.82 

-0.78 

-0.81 

-0.77 

-0.82 

0.63-0.69 

-0.91 

-0.86 

-0.80 

-0.73 

-0.81 

0.76-0.90 

0.93 

0.92 

0.76 

0.67 

0.79 

1.55-1,75 

-0.85 

-0.80 

-0.83 

-0,79 

-0.84 

2.08-2.35 

-0.,91 

-0.85 

-0.86 

-0.81 

-0,86 

Landsat MSS 

0. 5-0.6 

-0.82 

-0.79 

-0.81 

-0,76 

-0.81 

0.6-0. 7 

-0.90 

-0.85 

-0.81 

-0.74 

-0.82 

0. 7-0.8 

0.84 

0.84 

0.57 

0.46 

0.60 

0. 8-1.1 

0.91 

0.90 

0.77 

0.68 

0.79 
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Table A-2» The linear correlations of ratios and other transformations 
of reflectance measurements with agronomic characteristics 
of spring wheat canopies. 


Transformation 

Percent 

Soil 

Cover 

Leaf 

Area 

Index 

Fresh 

Biomass 

Dry 

Biomass 

Plant 

Water 

Content 

IR/Red* 

.89 

.90 

.76 

.67 

.78 

Red/IR 

-.94 

-.87 

-.80 

-.73 

-.81 

IR + Red 

.65 

.68 

.48 

.40 

.51 

IR - Red 

.95 

.93 

.81 

.72 

.83 

Vegetation Index 

.95 

.90 

.82 

•74 

.83 

Transformed Veg. Index 

.95 

.89 

.82 

.74 

.83 

Greenness 

.67 

.70 

.46 

. 36 

.48 

Brightness 

.90 

.89 

.68 

.58 

.70 


Red- 0.63-0.69 pm, IR - 0 .76-0.90 ym 
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variables. Similar results were reported by Tucker and Maxwell (1975). 

This region of low correlation reduces the ability of the 0.7~0.8 pm 
wavelength region to estimate crop canopy variables. The 0. 8-1.1 pm 
wavelength band has a greater correlation with each of the agronomic 
factors than the 0.7-0. 8 pm wavelength band, and the 0.76-0.90 pm 
wavelength band has a correlation with each agronomic factor that is 
greater then either the 0*7-0. 8 or the 0. 8-1.1 pm wavelength bands. 

Several investigators have examined the use of ratios, 
particularly infrared/red, and other transformations of two or more 
spectral bands to enhance the relationship between spectral and 
agronomic properties of vegetative canopies. Most of the investigations 
have been performed using Landsat MSS data, although the work of 
Tucker (1977) has been with field measurements of radiance. Compared 
to single wavelength bands these investigators have generally found 
stronger relationships using infrared/red rations, radiance differences, 
various vegetation indices, and the greenness-brightness transformation. 

Several transformations of the red and near infrared reflectances 
were evaluated in our study. These Included the ratios, differences, 
and sums of the two bands. In addition, transformations involving 
measurements in two visible and two infrared bands were evaluated. These 
included the vegetation index and transformed vegetation index developed 
by Rouse ^ al. , (1973) and the tasseled cap transformations into 
greenness and brightness developed by Kauth and Thomas, (1976). 

Several conclusions can be drawn from the results (Table A-2) of 
the analyses* The infrared/ red ratio, is sensitive to the amount of 
green Vegetation present, as is the red/lnfrared ration. The difference 
between the two bands is also highly correlated with the several 
measures of the amount of vegetation present, but their Sum is not. 

The correlations with the vegetation index and transformed vegetation 
index were similar to the ratios and differences. Geenness was also 
highly correlated with the canopy characteristics, but brightness was not. 

Compared to single bands, the use of the transformations, Infrared/red, 
red/infrared, vegetation index, transformed vegetazion index, and 



Figure A'- 6. Spectral reflectance of spring wheat canopies at several 

maturity stages. Measurements were made at Williston, North 
Dakota during May-August, 1976 and include plot? with 
different soil moisture levels, planting dates, nitrogen 
fertilization and cultivars. 




Maturity Stage 


Figure A- 7. Effect of maturity stage on the correlation of near 
Infrared reflectance (0.76-0.90 yin band) and 
agronomic characteristics of spring wheat canopies. 



greenness increased the correlation of the spectral data to th^ canopy 
variables, although not by a large amount, indicating that they are 
sensitive to the presence of green vegetation and may have the effect 
normalizing measurements' under varying conditions. However, the various 
transforms were all found to be highly correlated with each other and it 
is difficult to determine that any one is superior. 

Plant development and maturity (as opposed to growth or increase 
in size) cause many changes in canopy geometry, moisture content j and 
pigmentation of leaves which are also manifested in the reflectance 
characteristics of canopies. Figure A-6 shows the spectra of spring 
wheat at several different maturity stages. Maturity stage has a strong 
influence on the relationship of the canopy characteristics and 
reflectance as shown in Figure A-7. There is a definite decrease in 
the correlations during the heading-flowering and ripening stages. 
Learner et al. (1978) reported a similar decrease after heading in the 
correlation of spectral response with canopy variables for winter wheat. 

Prediction of Canopy Variables 

Estimation of canopy variables from multispectral data for use in 
crop growth and yield models is an important potential application of 
multispectral remote sensing. Understanding the relation of agronomic 
properties of crop canopies to reflectance in various regions of the 
spectrum leads to the development of regression models for estimating 
canopy variables from measurements in several wavelength bands. 

Table A~3 shows results using selections of one to six wavelength 
bands for prediction of canopy variables. By computing all possible 
regressions, the best subset of each size was selected considering 
the amount of variability explained and the bias of the resulting 
regression equation. The near and middle Infrared bands wore foiind to 
be most Important In explaining the variation In canopy variables. 

For leaf area Index and percent soil cover, the 0.76-0,90 pm wavelength 
band accounts for more of the variation than any other single band . 

The 2.08-2.35 pm wavelength band is the single most important band In 



Table A- 3. 


Selection of combinations of the best 1» 2, . . .6 wavelength 
bands for estimating percent soil cover, leaf area index, fresh 
biomass, dry biomass, and plant water content during the 
seedling thru flowering stages of crop development. 



No. 

Bands 



Bands Entered (um) 



Canopy 

r2 

* 0.45- 

0.52- 0.63- 0.76- 

1.55- 

2.08- 

Variable 

Entered 

Cp 0.52 

0.60 0.69 0.90 

1.75 

2.35 


Percent 

1 

.86 

132 




X 



Soil 

2 

.92 

16 




X 


X 

Cover 

3 

.92 

15 

X 



X 


X 


4 

.93 

4 

X 

X 


X 


X 


5 

.93 

5 

X 

X 

X 

X 


X 


6 

.93 

7 

X 

X 

X 

X 

X 

X 

Leaf 

1 

.84 

37 




X 



Area 

2 

.87 

7 




X 

X 


Index 

3 

,88 

2 




X 

X 

X 


4 

.88 

4 

X 



X 

X 

X 


5 

.88 

5 

X 

X 


X 

X 

X 


6 

.88 

7 

X 

X 

X 

X 

X 

X 

Fresh 

1 

.73 

239 






X 

Biomass 

2 

.76 

211 


X 




X 


3 

.83 

109 


X 

X 



X 


4 

.88 

41 


X 

X 

X 


X 


5 

.90 

12 

X 

X 

X 

X 


X 


6 

.93 

7 

X 

X 

X 

X 

X 

X 

Dry 

1 

.65 

252 






X 

Biomass 

2 

.67 

229 

X 

X. 






3 

.81 

78 


X 

• X 



X 


4 

.84 

44 


X 

X 

X 


X 


5 

.87 

20 

X 

X 

X 

X 


X 


6 

.88 

7 

X 

X 

X 

X 

X 

X 

Plant 

1 

.75 

201 



t 



X 

Water 

2 

.77 

175 


X 




X 

Content 

3 

.83 

98 

X 

X 


X 




4 

.88 

34 


X 

X 

X 


X 


5 

.90 

9 

X 

X 

X 

X 


X 


6 

.90 

7 

X 

X 

X 

X 

X 

X 


*The regression equation is unbiased when the ’Cp’ value is equal to or 
less than the number of terms in the equation. 


explaining the variation in fresh biomass, dry biomass, and plant water. 
The 2. 08-2.; 35 pm wavelength band is one of the two most important 
bands in explaining the variation in percent soil cover and one of the 
three most important bands explaining the variation in leaf area index. 

From Table A-3 the difference between the number of bands entered 

2 

that produce a near maximum R and the number of bands entered where 

the resulting prediction equation is unbiased can also be examined. 

An unbiased equation results when the value is equal to or less 

than the number of terms in the resulting regression equation (Mallows, 

1973). For leaf area index and percent soril cover, the near maximum 
2 

R value is reached after the entry of three out of the six possible 
thematic mapper bands. However, the values indicate that five 

bands would have to be used to have an unbiased prediction of leaf area 
index and four bands would be necessary for percent soil cover. 

The agreement between the measured and predicted percent soil 
cover, leaf area index, and fresh biomass is shown in Figure A-8. 

Similar results were obtained for the other canopy variables. There 
are several factors that make a perfect prediction impossible for these 
data, including: (1) the agronomic measurements of the crop canopy 

were subject to measurement error, (2) plant maturity stage has an 
effect on reflectance (for example, a canopy with an LAI of 1.0 early in 
the season has a different spectral response than a canopy with a 
similar LAI later in the season) , (3) the data that the prediction 
equations are derived from contains variation induced by the different 
agronomic treatment levels , and (4) the time of day that the data were 
collected may have some effect on canopy reflectance.’ Despite the 
variation induced by each of these factors, measurements in a small 
number of wavelength bands in important regions of the spectrum can 
explain much of the variation in canopy variables and thus, result in 
satisfactory predictions of canopy variable-s . 

2 

Table A~4 shoves the maximum R value obtainable using the Landsat 
bands, the best four out of the six posslbJe thematic mapper bands, and 
then all six thematic mapper bands to predict each canopy variable. In 
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Figure A-8. Comparison of measured and predicted percent soil 

cover, leaf area index, and fresh biomass of spring 
wheat canopies. 
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Table A-4. The R values for predictions of percent soil cover, leaf 
area index, fresh and dry biomass, and plant water content 
with four Landsat MSS bands, the best four thematic mapper 
bands, and the six thematic mapper bands. 


Wavelength 

Bands 

Percent 

Soil 

Cover 

Leaf 

Area 

Index 

Fresh 

Biomass 

Dry 

Biomass 

Plant 

Water 

Content 

Landsat MSS Bands 

0.91 

0.86 

0.86 

0.84 

0.85 

Best Four Thematic 

Mapper Bands 

0.93 

0.88 

0.88 

0.84 

0.88 

Six Thematic Mapper 

Bands 

0.93 

0.88 

0.91 

0.88 

0.90 


See Table A- 3. 


every case die best four out of six tliematic mapper bauds explained 

more of the variation in a canopy variable tlien the four l.andsat MSS 

bands. Addition of the other two thematic mapper bands resulted in only 

2 

small increases in the R*" values. 

1.3 Summary and Goneluslons 

Spectral and agronomic measurements of spring wheat canopies were 
analyzed to determine the relation of agronomic' properties of crop 
canopies to their spectral reflectance. A strong relationship between 
spectral response and percent soil cover, leaf area index, biomass and 
plant water content was found. The relationship, however, is influenced 
by crop maturity. The best time period for assessing these canopy 
vax'iables is from the tillering to heading stages of development. Prior 
to tillering the spectral response is strongly dominated by the soil 
background and, as the crop begins to ripen, the spectral sensitivity 
to measures such as leaf area index, biomass, and plant water content 
decreases. 

In each wavelength region, the correlation of the thematic mapper 
band with crop canopy variables was greater than that of the corresponding 
Landsat MSS band. Regression equations developed to explain the variation 
in crop canopy variables showed that the 2.08-2.35 pin wavelength band 
was the single most important band in explaining the variation in fresh 
biomass, dry biomass, and plant water content; whereas, the near infrared 
band (0.76-0.90 pm) explained the most variation in leaf area index and 
percent soil cover. The results demonstrate the importance of collecting 
spectral Information in the middle infrared wavelength region, as well 
as the visible and near infrared, for crop assessments. 

.2 

The R values for comparisons of measured and predicted canopy 
variables ranged from 0.80 to 0.91 when three or more spectral bands were 
included, indicating the potential for using remotely sensed spectral 
measurements to characterize the status of crops. Analyses showed that 
the best four thematic mapper bands could estimate crop canopy variables 
better than the four Landsat bands. The difference is attributed to the 
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narrower and more optimum placement of the thematic mapper bands in 
relation to the spectral characteristics of vegetation. 


The strong relationship between spectral reflectance and different 
crop canopy variables illustrates the potential to monitor crop growth 
and predict yield. Future research needs to investigate the amount of 
variation iisJuced by different agronomic treatment factors on spectral 
reflectance and whether important treatment factors are spectrally 
separable. The type of prediction equations developed in this 
investigation need to be extended to several years of data, then used 
to estimate independent data sets. 
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2. EFFECTS OF CULTURAL AND ENVIRONl-IENTAL FACTORS ON THE REFLECTANCE 
CHARACTERISTICS OF SPRING WHEAT CANOPIES 

The crop canopy is a dynamic entity influenced by many cultural 
and environmental factors. Therefore, it is important to quantify and 
understand the sources of variation in spectral measurements of crops. 

Some of the variation may be associated with important agronomic factors 
which it may be desirable to inventory or monitor (for example, dryland 
vs. irrigated wheat). On the other hand it is also important to know 
the magnitude of variation associated with a factor such as cultivar 
which We would most likely not want to identify or monitor. 

Examples of spectra acquired in 1976 at the Williston, North Dakota, 
Agriculture Experiment Station are shown in Figure A-9 to illustrate some 
of the effects of agronomic treatments on the spectral response of spring 
wheat. Further experiments using an improved experimental design were 
conducted in 1977 to determine the effects of the various agronomic 
treatments (soil moisture availability, planting date, nitrogen fertiliza- 
tion, and cultivar) on the reflectance of spring wheat canopies. The 
results of the 1977 investigation are reported in this section. 

2. 1 Experimental Approach 

Data were collected at the Williston, North Dakota Agriculture 
Experiment Station during the summer of 1977. The Exneriment Station is located 
in the gently rolling uplands above the Missouri River Valley and is 
representative of dryland farms of the region. The average annual precipi- 
tation of 36 cm is just adequate for successful farming operations in 
normal years. The area is semi-arid and humidity remains low throughout 
the year. 

The soil in the study area is a Williams loam (Typic Argiustoll) 
which was developed in glacial till. The A^^ surface horizon is dark 
brown (lOYR 3/2) when moist, but very light (lOYR 4/3) in color when dry. 

The predominant crop grown in the northwest area of the state is hard red 


Bidirectional JRefi&clance Factor [%) 



Wavelength (ym) 


Figure A-9. Effects of agronomic treatments on the spectral reflectance 
of spring wheat. Spectra were measured on June IS, 

1976, during the stem extension stage of development, 
except for the spectra of cultivars which were measured on 
July 16 after heading. 
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spring wheat. Because of the limited amount of rainfall, the majority 
of land is planted to crops every other year and is left fallow in 
Intermediate years to allow subsoil moisture to accumulate. 

The spring wheat experiment under investigation was a split plot 
design. Within each available soil moisture condition (whole plot), 
there were two replications of a factorial design with cultivar, nitrogen 
fertilization, and planting date as experimental treatments! 

Available Soil Moisture 

Low (wheat crop in 1976) 

High (fallow in 1976) 

Planting Date 

Early (May 9, 1977) 

Late (May 23, 1977) 

Cultivar 

Semi-dwarf (Olaf) 

Standard (Waldron) 

Nitrogen Fertilization 
None 
44 Icg/ha 

Spectral measurements were made with an Exotech Model 100 Radiometer 
in four wavelength bands (0.5-0. 6, 0.6-0. 7, 0.7-0. 8, and 0.8-1. 1 pm). 
Measurements in all bands were taken instantaneously and recorded by a 
hard copy data logger, A boom mounted on a van supported the radiometer 
at 3.5 m above the canopies and 3.5 m away from the van. At this 
elevation and with a 15° field of view the sensor viewed a 0.9 m diameter 
ground area. The reflectance of all 32 plots was measured at approximately 
hourly intervals during the middle part of the day when the sun angle 
was greater than 45 . 

A 1.2 m square standard painted with highly reflecting barium 
sulfate was used as a basis for the bidirectional reflectance factor 
calibration. The reflectance standard provided a field calibration 
reference with stable, known reflectance properties. A dark level 
response of the instrument was also obtained by holding an opaque, light- 
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tight appgiratus against the instrument's entrance ports to measure the 
internal system noise or deviation from ?ero. The response of the 
calibration panel was measured approximately every 15 minutes during data 
collection over the plots and the dark level every 30 minutes. These 
values were then used in the following equation to calibrate readings 
taken over the plots; 


= DrW : draj • 

Where, BRF(X) - bidirectional reflectance factor (%) at a specific wavelength 
interval (A), 


Ds(X) = response of Instrument to scene (plot), 

ds(X) = dark, level response of instrument taken closest in time to 
scene, 

Dr(A) « response of instrument to painted barium sulfate reference 
standard. 


dr(X) » dark level response of instrument; taken closest in time to 
reference standard measurement, and 

Rr(X) * reflectance of painted barium sulfate reference standard 

(measurement in laboratory by comparison with pressed barium 
sulfate) . 


The plots were planted in 17.8 cm (seven inch) rows and were 
rectangular in shape (3. 5 x 15.3 m) . Leaf area index, fresh and dry 
biomass, growth stage, plant height, number of leaves per plant, percent 
of green, yellow and brown leaves were measured on a one meter row sample 
from each plot on each day that reflectance measurements were made. 

In addition to the daily meteorological records kept by the Experiment 
Station, a compact weather station was set up each day that spectral 
data were collected. Meteorological measurements including air temperature, 
barometric pressure, relative humidity, wind speed and direction, and 
total irradiance were recorded continuously through the day on strip charts. 

The objective of the data analysis was to determine the significance 
of the agronomic treatments for each of the dates (crop maturity stages) 
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and wavelength bands. Since previous work (Crecelius, 1978) has 

indicated that there is some diurnal variation in the spectral reflectance 

of wheat canopies, a covariance analysis was used to adjust the tfeatment 

means of the reflectance variables to be estimates of what they would 

be if all plots were measured at the same time (Steel and Torrie, 1960) . 

A separate covariance analysis was run for each date and wavelength band 

2 

with time and time as covariates to adjust for the two components of 
diurnal variation which are synraietric and asymmetric about solar 
noon (Crecelius, 1978), 

2.2 Results and Discussion 

Initial analysis of the data indicated that the whole plot error 
(between soil moisture levels) was consistently larger than the split 
plot error for each date and hand analyzed. It was determined that the 
replications in this experiment were not a result of random selection 
but described a location within the whole plot. Thus, differences due 
to available soil moisture could not be tested statistically. 

A summary of the significant results based on the covariance analysis 
is presented in Table A-5 for each date and wavelength band. The effects 
indicated there had a significant (a=, 01) effect on spectral response 
after adjusting for the symmetric and asymmetric components of time. 

Planting date was generally significant for all dates and at all wave- 
lengths. Although no statistical test for the main effect of soil 
moisture level was available, significant interactions of other factors 
with soil moisture provided indications of its significance. An additional 
method of ascertaining the overall importance of the experimental treatments 
is to consider the percentage of total variability accounted for by each 
treatment (Table A-5). This approach permits an evaluation of the effect 
of available soil moisture. 

On June 1, planting date is the primary factor affecting spectral 
response (Table A-5). With a 14 day difference between the early and 
late planting dates there are differences in the agronomic, as x<rell as 
spectral, characteristics of the wheat canopies (Table A-6) . The wheat 
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Table A-5. Percent of total variation In reflectances , greenness and 
brightness accounted for by soil moisture level, and 
significant agronomic factors and interactions. 


Mission 

Date 

Treatment , 
Effect^ 

Wavelength Band (pm) 

Transformations 

.5-. 6 

.6-. 7 

.7-. 8 

.8-1.1 

Greenness Brightness 

June 1 

M 

4. 7 

3.6 

9.4 

8.3 

4.6 

8.2 


P 

58.6 

34.4 

81.9 

83.5 

42.1 

79.2 


MP 

5.4 

8.4 

- 

- 

10.1 

0.5 


CNP 

- 


0.3 

- 

- 

- 


MNP 

- 

- 

0.3 

- 

- 



MCP 

- 

- 

0.2 

- 

- 

- 


MCN 

0.9 

- 

— 

- 

- 

- 

June 18 

M 

5.7 

5.9 

0.2 

0.9 

3.9 

0.1 


N 

2.5 

- 

i - 

— 

■ — . 

.. 


P 

19.3 

22.7 

j 26.9 

41.5 

47.9 

14.6 


MN 


■- 

1 3.2 

3.6 


2.2 


MP 

39.0 

32.7 

3.2 

- 

11.6 

10.1 


MNP 

3.3 

3.3 

- 

- 

- 

- 

July 14 

M 

52.4 

60.5 

9.2 

53.4 

61.4 

0.0 


q 

1.2 

1.0 

- 

3.0 

2.1 

- 


N 

1.9 

1.6 

5.5 

5.5 

3.8 

- 


P 

6.5 

4.1 

- 

4.7 

6.1 

- 


MP 

5.3 

6.0 

12.4 

- 

0.8 

21.1 


MCNP 

— 


- 

2.2 

1.1 

- 

July 20 

M 

64.6 

51.9 

6.9 

34.7 

62.8 

15.0 


C 

- 

- 

- , 

- 

2.8 

- 


N 

- 

- 

- 

- 

3.1 



P 


26.7 

- 

- 

21.0 

- 


MP 

5.5 

■ 

25.2 

- 

- 


August 5 

M 

61.8 

53.4 

33.9 

4.7 

58.2 

40.8 


C 

5.7 

12.1 

- 

- 

3.4 

- 


N 

- 

- 

- 

- 

3,1 

- 


P 

5.7 

12.1 

- 

- 

12.3 

5.3 


MP 

0 . 3 


" 

7.5 




Total variation Is considered that due to treatments, diurnal Variation, 
and experimental error. ’ 


M, available soil moisture; P, planting date; N, nitrogen fertilization; 
C, cultivar. 


Table A-6. Mean agronomic characteristlrq of cn-r':.,r. w ^ 

level and planting date (Williston, North cSotarSyy) moisture 


Soil 

Measurement Jfcisture Planting Ifatnrlty Plant Soil Aral 
Dai_e Stage Height Cover Index 


June 1 

High 

Early 

2.1 

14 

7 

0.7 



Late 

1.0 

Q 

2 

0.1 


Low 

Early 

2.1 

13 

6 

0.3 



Late 

1.0 

9 

2 

0.2 

June 18 

High 

Early 

2.4 

29 

47 

1 7 



Late 

2.2 

23 

23 

0.7 


Low 

Early 

2.4 

24 

20 

0.8 



Late 

2.2 

23 

22 

0.7 

July 14 

High 

Early 

4.7 

64 

52 

1 3 



Late 

4.4 

58 

60 

2.1 


Low 

Early 

4.7 

46 

42 

0 5 



Late 

4.6 

46 

49 

0.8 

July 20 

High 

Early 

4.8 

67 

50 

0 7 



Late 

4.7 

68 

51 

1.2 


Low 

Early 

5.1 

45 

38 

0 2 



Late 

4.8 

48 

44 

0.5 

August 5 

High 

Early 

5.3 

65 

46 

0 0 



Late 

5.2 

66 

46 

0.1 


Low 

Early 

5.3 

46 

30 

0 0 

■ ' 


Late 

5.2 

45 

35 

0.0 


Fresh 

Biomass 

Dry 

Biomass 

Plant 

Water 

Content 

g/m^ 

g/m^ 

g/m^ 

72 

8 

64 

32 

1 

31 

50 

7 

43 

37 

3 

34 

769 

105 

664 

335 

30 

306 

407 

50 

358 

315 

29 

285 

1528 

620 

908 

1709 

480 

1229 

723 

2S5 

438 

715 

215 

500 

1525 

603 

922 

1454 

476 

978 

781 

307 

480 

694 

224 

470 

672 

466 

206 

1117 

550 

567 

439 

237 

202 

567 

259 

308 
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planted at the later planting date (May 23) is at the seedling stage of 
maturity, while the early planted wheat (May 9 planting date) Is beginning 
to tiller. The differences In spectral response are attributed to the 
differences In percent soli cover, leaf area index, and biomass, flantlng 
date continues to be the primary agronomic factor affecting reflectance on 
June 18 during the later phases of tillering. 

By July 14 differences in vegetative growth and spectral response 
due to planting date are reduced as the later planted wheat, growing at 
a faster rate due to warmer temperatures and longer days during its 
early stages of development than the early planted wheat, is nearly the 
same size as the early planted wheat. At this time, heading to flowering 
stages of maturity, soil moisture availability is the major factor 
affecting the growth and development (Table A-6) and spectral response 
(Ta1>le A*^5) of the wheat canopies. Although there are significant 
differences due to culLlvar and nitrogen fertilization, they are small 
compared to those associated with available soil moisture. 

A week later, on July 20, soil moisture level is the most important 
factor accounting for differences in spectral response (Table A-5). 

Effects due to cultlvar and fertilization are not significant except in 
greenness. 

On August 5, as the canopies are ripening, soil moisture level 
continues to be the primary factor influencing the spectral response. 

Since the canopies have similar maturity stages, the differences in 
spectral response are attributed to the differences In percent soil cover 
and biomass (Table A-6). 

2.3 Summary and Conclusions 

The effects of four agronomic treatments (soil moisture, planting 
date, nitrogen fertilization, and cultlvar) on the reflectance of spring 
wheat canopies was investigated. Early in the growing season planting date 
is the primary agronomic factor influencing the reflectance of spring wheat 
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The spectral differences are primarily due to differences in the amount of 
vegetation present. Later in the season at the heading to ripening stages, 
the level of soil moisture becomes the .most important factor. Wheat grown 
on land with higher levels of available soil moisture had a greater 
percent soil cover, leaf area index, and biomass causing increased near 
infrared and reduced, visible reflectance. In these experiments cultivar 
and nitrogen fertilization had relatively little effect on the spectral 
response of spring wheat. Again, examination of the agronomic data 
shows that these two treatments had little influence on the growth and 
development of the spring wheat canopies. The primary difference in the 
two cultivars was in plant height, rather than in leaf area or biomass. 

And, since the soils of this area of North Dakota are relatively high in 
nitrogen supplying capacity, the addition of nitrogen fertilizer had only 
minor effects of the growth of wheat. 
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3. ASSESSMENT OF CROP DISCRIMINABILITY AS A FUNCTION OF CROP 
MATURITY STAGE AND WAVEI.ENGTH BAND SELECTION 

The application of remote sensing to crop inventories relies upon 
the ability to detect and identify the crops of interest from the spectral 
data. Unless correct crop identification can be consistently made, the 
goal of accurately predicting crop acreage and production cannot be achieved. 

The overall objective of these analyses was to assess the potential 
for crop discrimination. The specific objectives were; (1) to 
determine which cover types may be confused, (2) to determine at what 
growth stages separability is maximized, and (3) to determine which regions 
of the spectrum are best for crop discrlminabillty . 

3.1 Experimental Approach 

LACIE Field Measurements data acquired by the helicopter spectrometer 
over the intensive test sites at Finney County, Kansas, In 1975 and 
Williams County, North Dakota, in 1976 were analyzed. These data are xirell- 
suited for this Investigation because of the frequent measurements made in 
narrow wavelength bands and the associated ground truth observations. 

Spectral data were analyzed, using the Landsat MSS and proposed thematic 
mapper wavelength bands, as well as the greenness and brightness transfor- 
mations of reflectance. 

Spectral plots of the data were made in two ways: (1) several crops 

on the same measurement date and (2) a single crop over several dates. 

The plots gave qualitative indications of what crops may be separable at 
any given time and also show changes in a crop over the growing season. 

An extensive set of example plots can be found in the report entitled, 

"Crop Spectra from LACIE Field Measurements," (Hixson et al. , 1978). 


Analysis of variance and range tests were used to determine dates and 
wavelength bands where significant differences in spectral response occurred. 
Unitemporal and multitemporal multivariate discriminant analyses were 
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performed to investigate crop discriminability . Perceitt correct 
classification and percent commission errors were computed. 

3.2 Results and Discussion 

Williams County, North Dakota 

In general, 98 percent of v;heat spectra were correctly classified in 
the June 26 and later data for both the Landsat MSS and thematic mapper 
bands (Table A- 7 ) . On the first three dates, the percent of wheat spectra 
correctly identified was generally above 80. Spectra of fallow fields 
were almost equally well identified on most dates with both Landsat MSS and 
thematic mapper bands. Fifty percent of the spectra of pasture were confused 
with fallow after July 20 when using Landsat MSS bands, but substantially 
less confusion occurred when using the thematic mapper bands. 

There was Some variation in the percent of wheat spectra correctly 
classified based on the size of the training set. The results presented in 
Table A-7 are based on using all available spectra for training. Using 
a random sample of either 50 or 25 percent of spectra produced very similar 
results , almost identical on June 25 and later dates. 

Although the amount of training used did not significantly alter the 
percent correct classification achieved, the training method used did 
have some effect. Figure A-IO presents results obtnincd from sampling 
about 50 percent of the spectra using sampling of fields rather than sampling 
of Individual spectra. The percent of wheat correctly identified, particularly 
in the Landsat MSS bands, was decreased from that obtained using 50 percent 
of the spectra randora3.y sampled. 

When cumulative spectral information from several dates was used, 
higher classification accuracies were achieved. The percentage of wheat 
spectra correctly classified ranged from 97 to 100 for both Landsat and 
thematic mapper bands. The confusion of pasture with other cover types 
was not enCounteted for either Landsat MSS or thematic mapper bands. 
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Table A-7. Percent correct classification of spring wheat, fallow, 
and pasture using unitemporal information (Williams Co. 
ITS, 1976). 


Mission 

Date 

Land. 

sat MSS Bands 

Thematic Mapper Bands 

Wheat 

Pasture 

Fallow 

Wheat 

Pasture 

Fallow 

May 13 

96.3 

90.6 

96.7 

95.7 

96.6 

96.7 

May 28 

76.4 

81.7 

64.4 

84.2 

85.9 

72.8 

J une 1 7 

90.5 

60.0 

93.9 

96.4 

83,3 

91.7 

June 25 

98.8 

73.1 

96.2 

99.6 

82.8 

95.7 

July 6 

98.1 

84.0 

99.4 

98.9 

93.3 

99.1 

July 20 

99.6 

48.1 

98.9 

99.6 

88.0 

98.9 

July 28 

99.2 

38.1 

97.0 

98.8 

35.8 

98.3 

August 9 

97.0 

96.4 

99.0 

98.0 

100.0 

99.2 

August 19 

98.0 

57.6 

96.6 

98.4 

98.0 

98.8 


Table A-8. Percent correct classification of spring wheat, fallow, 

and pasture using the Landsat MSS bands and multitemporal 
Information (Williams Co. ITS, 1976). 


Combination of 


Mission Dates 

Wheat 

Pasture 

Fallow 

6/17, 7/6,-7728 

99.3 ‘ 

92.2 

99,4 

5/28,6/25,7/28 

98.8 

88.0 

98.8 

5/28, 6/25, 7/20, 8/9 

99.5 

100.0 

99.4 

6/17, 6/25, 7/6, 7/28 

100.0 

95.6 

100.0 
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Figure A-10. Comparison of Landsat MSS and thematic mapper spectral bands for spring wheat 
identification (Williams County, North Dakota, 1976). 
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An investigation was then carried out to find a reduced set of 
dates which, when combined, would give results comparable to using the 
full data set. The combinations of dates used were selected to include 
Information from several important grov/th stages. Again, using raultitemporal 
information, the confusion of pasture with other cover types was not encountered 
and all cover types were identified with high accuracy (Table A-8 ). 

Similar uniiteraporal analyses using the brightness and greenness 
components of reflectance were also conducted. The results of these 
analy,ses did not show an improvement in classification of wheat spectra 
over the use of two Landsat MSS bands. 

Finney County, Kansas 

Univariate analyses of variance showed that for every wavelength 
band and for every date except (.7-. 8 ym on June 17), there was a signi- 
ficant difference in reflectance due to crop type. All data were used for 
training and evaluation of the classification. The results are shown in 
Table A- 9 . 

It appears that wheat and alfalfa can be well distinguished in 
November. Some confusion exists between these cover types and fallow, 
perhaps due to the small proportion of soil cover or insufficient statistics 
to define the spectral distribution of fallow. 


In March and April, most wheat has reached a stage of greenness 
where it can be differentiated from fallow and these cover types remain 
separable on all following measurement dates. A few wheat fields are 
still being missed, but almost no fields are being falsely identified as 
wheat. 

On May 14, alfalfa and wheat are distinctly differentiable and remain 
so on all following measurement dates. Corn, however, has not yet reached 
a growth stage where it is separable from fallow land. Com and fallow 
become spectrally separable by May 21 and then remain so. Beginning on 
June 8, sorghum can be correctly identified with the exception of some 
fallow land being classified as sorghum. Fallow, which consists of bare soil 


Table A-9. Summary of the separability of cover types at several times 
during the 1974-75 growing season for the Finney County, 
Kansas, intensive test site. 


Overall Percent 
Correct Identification 


Mission 

Date 

Classes 

Landsat 

MSS 

Thematic 

Mapper 

Comments 

Nov. 5 

wheat, alfalfa, 
fallow, other 

77 

76 

Some wheat was classi- 
as alfalfa and fallow 

March 20 

wheat, alfalfa, fallow 

87 

88 

Some wheat was missed. 

April 8 

wheat, alfalfa, fallow 

88 

87 


May 14 

wheat, corn, fallow, 
alfalfa, other 

77 

76 

Major problem is 
being classified as 
another cover type. 

May 21 

wheat, corn, fallow, 

89 

90 

Same as May 14. 

June 2 

wheat, corn, fallow, 
alfalfa, sorghum, 
other 

75 

79 

Wheat and alfalfa 
correctly identified. 
All other cover types 
are confused. 

June 9 

wheat, corn, fallow, 
alfalfa, sorghum, 
other 

67 

88 

Wheat and sorghum 
correctly identified by 
the Landsat bands ; 
alfalfa was also 
correct with thematic 
mapper. 

J une 1 7 

wheat, corn, fallow, 
alfalfa, sorghum, other 

78 

83 

All crops are separable 
from fallow; wheat and 
alfalfa separable; 
other crops confused. 

J une 26 

wheat, corn, fallow. 

88 

90 

Same as June 17. 

July 6 

wheat, com, fallow, 
alfalfa, sorghum, other 

84 

92 

Same as June 17 except 
more confusion of 


sorghum and fallow. 
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and several degrees of stubble, is confused with other crop types early 
in the season. 

On June 9 and July 6, there is a noticeable difference in crop discri- 
minability between performance of the Landsat MSS and proposed thematic 
mapper bands with accuracies of thematic mapper bands being higher. 

3.3 Summary and Conclusions 

Discriminant analyses of crop discri^nability were conducted using 
FSS data with simulated Landsat MSS and thematic mapper wavelength bands. 

The analyses included data from the Williams County, North Dakota, and 
Finney County, Kansas, intensive test sites. 

Spring wheat in North Dakota was confused with pasture by the Landsat 
MSS bands when only single date Information was used. If raultitemporal 
information was used, however, the two cover types were separable. Training 
method was found to have a greater effect on classification results than 
amount of training. Insufficient data were acquired at the test site for 
quantitative analyses describing the separability of wheat from other 
small grains. 

Winter wheat in Kansas was occasionally confused with alfalfa and 
was confused with fallow land early in the season before sufficient soil 
cover had been reached for the wheat to appear characteristically green. 

Corn and sorghum could not be identified until much later in the season due 
to their late planting date. 

Classification performance for the thematic mapper bands was somewhat 
higher than for the Landsat MSS bands. Use of the brightness and greenness 
components of reflectance did not show an Improvernent in classification of 
wheat spectra, over the use of two Landsat bands. 

3.4 References 
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B. Field Measurements Data Management 


The development of the field research data library at Purdue/LARS 
was initiated in tihe fall of 1974 by NASA/ Johnson Space Center (JSC) 
with the cooperation of the United States Department of Agriculture (USDA) 
as a part of the Large Area Crop Inventory Experiment (LACIE) . The 
purpose for developing the data base is to provide fully annotated 
and calibrated multitemporal sets of spectral, agronomic, and meteorological 
data for agricultural remote sensing researchers. Spectral, agronomic, 
and meteorological measurements oVer primarily wheat have been made 
on three LACIE test sites in Kansas, North Dakota, and South Dakota 
for three years. During this past year the data library was expanded to 
include data collected for com and soybean experiments in Indiana. 

Milestones reached during this year have been: (1) completion of 

the 1976-77 data processing, (2) development of graphical and 
statistical analysis software, (3) development of documented data 
processing software, and (4) distribution of data to researchers represent- 
ing university, government and commercial research organizations. 

A report entitled "Crop Spectra from LACIE Field Measurements" 
was prepared and distributed this part year. This document contains 
examples of the spectrometer data in the form of spectral reflectance 
curves Illustrating the major sources of spectral variation of wheat and 
related crops. The examples include variations in wheat spectra due 
to differences in maturity, biomass, soil color, and soil moisture, as 
well as comparisons of spectra of wheat and other crops. 


' ' 

This section describing the results of Task IB,- Field Measurements Data 
Management, was prepared by the task leader, Larry Biehl. The contributions 
of Jeanne Etheridge, Cathy Kozlowski, Jerry Majkowski, and Jeff McMeekin to 
the data processing and analysis software development and Cathy Axtell and 
staff of the LARS Computer Facility for data preparation are gratefully 
acknowledged. 


1. FIELD RESEARCH DATA LIBRARY AND DISTRIBUTION 


The general organization of the field research data library is 
illustrated in Figure B-1. The data in the library includes spectral 
measurements (multispectral scanner and spectrometer), agronomic 
measurements, meteorological measurements, photography, mission logs 
and data verification reports. The data formats available to researchers 
are digital tape, film, and data listings. 

« 

The Instruments used to collect spectral data are of two types, 
multispectral scanners and spectrometers/radiometers . The multispectral 
scanners include the 4-channel Landsat MSS, the 11-channel MMS on the 
NASA-T'3 aircraft, and the 24-channel MSS on the NASA-C130 aircraft. 

The spectrometers include the NASA/ERL Exotech 20D field system, the 
NASA/JSC Field Spectra Acquisition System (FSAS) , the NASA/JSC Field 
Spectrometer System (FSS) , the Purdue/LARS Exotech 20C field system. 

The radiometer includes the Purdue/LARS Exotech 100 field system. 

The spectrometer data are processed into a comparable format 
(spectral bidirectional reflectance factor) in order to make meaningful 
comparisons of the data acquired by the different sensors at different 
times and locations. The spectrometer tapes contain the spectral 
bidirectional reflectance factor measurements along with most of the 
corresponding agronomic and meteorological measurements. The data is 
stored on the tape in LARS spectrometer data storage tape format. 

The multispectral scanner data stored on the LARSYS Version 3 
formatted tapes are approximately linearly related to scene radiance. 

The information is available for the researcher to calibrate the 
scanner data to in-band bidirectional reflectance factor if he wants 
to. Some of the multispectral scanner data is available in universal 
format. 

In the past twelve months, 45 aircraft scanner runs and 30,000 
spectrometer runs have been processed and made available in the field 



Figure B-1. Organization of Field Research Data Library. EXOSYS and LARSYS are Purdue/LARS software 
systems to analyze spectrometer and multispectral scanner data. 
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research data library. The status of data processing is sununarized in 
Tables B-1 and B-2. Processing of the remaining 1976-77 crop year 
data, along with the sun angle-view angle-view direction data, was 
completed earlier this year. Processing of the 1977-78 data was 
begun this fall. 

The Field Research Data Library Catalog summarizes the data 
available. The catalog is divided Into separate volumes - one for 
each crop year during which data were collected. Volume III, 1976-77, 
of the catalog was updated this past year. Discussions with the past 
and present users of the catalog and data have provided direction for 
future updates of the catalog. 

Six institutions, representing university, government, and 
commercial research organizations, have requested and received the 
LACIE field research data this year. In addition all the data are 
routinely available to researchers at Purdue/LARS and investigators 
at NASA/Johnson Space Center have direct access to the digital data 
via the remote terminal to the LARS' computer located there. 

Another organization Inquired about specialized processing of 
some of the aircraft scanner data into thematic mapper resolution, 
but did not formally request it. Table B-3 summarizes the data 
distribution during the past year. 


Table B-1. Summary of LACIE Field Measurements data in the 


Field Research Library by year 
data type. 

, instrument* 

and 



Crop Year 

-• 

Instrument /Data Type 

1975 

1976 

1977 

Landsat MSS 



.1 ... ^ 

Whole Frame CCT (Frames) 

20 

62 

34 

Aircraft Scanner 
(Dates /Runs) 

19/149 

16/97 

7/35 

Helicopter Mounted Field 
Spectrometer (Dates/Runs) 




Field Averages 

19/2,343 

27/2,193 

21/1873 

Individual Scans 

19,29,579 

27/38,476 

21/34.060 

Truck Mounted Field 
Spectrometer (Dates/Runs) 




FSAS 

6/65 

23/322 

16/426 

Exotech 20C 

24/1599 

18/2542 

21/1075 

Exotech 20D 

45/645 

- - _ — 

— — — 
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Table B-2. Status of Preparation of 1977-78 Multicrop Supporting 
Field Research Data. 


Sensor System/Data Type 

Completed & 
In Library 

In 

Processing 

At 

JSC 


Dates /Observations 

Landsat MSS 

Whole Frame CCTS 

2 

0 


Aircraft Scanner 

3/15 

1/5 

- - - 

Helicopter-Mounted Field Spectrometer 

Field Average Data 

0 

6/ 

2/ 

Individual Scans 

0 

6/ 

2/ 

Truck Mounted Field Spectrometer 

2/58 

22/1097 

- - - 


The 1978 Landsat tapes for the Hand County site are in processing at GSFC. 
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Table B-3. Recipients of Field Research data during 1977-78. 


Institution 

Data Requested and Received 

Enviromental Research 
Institute of Michigan (ERIM) 

Spectrometer data tapes, ground 
photos, flight logs, field maps. 

Texas A&M University, 
Remote Sensing Center 

Spectrometer data tapes, ground photos 
flight logs, field maps 

NASA, Goddard Space Flight 
Center 

Spectrometer data tapes, Landsat 
MSS data tapes, aerial photography, 
ground photos, field maps. 

General Electric Corporation 
(contract with NASA/GSFC) 

Spectrometer data tapes, agronomic and 
meteorological measurements, 
field maps, ground photos, aerial 
photography. 

General Electric Corporation 

Landsat MSS data tapes, agronomic 
measurements, Landsat color 
composite imagery. 

USDA, Agriculture Research 
Service, Weslaco, Texas 

Helicopter (FSS) photography, 
agronomic measurements, 
meteorological measurements, 
field maps. 

University of South Florida, 
Department of Mathematics 

Spectrometer data tapes. 
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DEVELOPMENT OF GRAPHICS, STATISTICS AND DATA PROCESSING SOFTWARE 


During the LACIE Field Measurements Project over 100,000 
observations of spectrometer data were collected. Additional observa- 
tions were collected this year in preparation for multicrop studies. 

These increasingly large numbers of field measurements type observa- 
tions have brought about a need for more powerful research tools to 
verify and analyze the data efficiently and effectively. During this 
past year, several new analysis software tools were developed and 
Implemented on the Purdue/LARS IBM 370/148 computer. These include 
new graphics capabilities and additional statistical capabilities 
(such as clustering) for spectrometer data. Also, an updated version of 
SPSS (Statistical Package for the Social Sciences) having additional 
statistical analysis capabilities was obtained and made available to 
researchers on the Purdue/LARS computer. 

The software system implemented on the Purdue/LARS computer to 
access the spectrometer data and associated agronomic and meteorological 
measurements stored on tape is called EXOSYS. The version which 
contains the new graphics and statistics capabilities is termed EXOSYSDV. 
The new analysis tools for researchers include capabilities to: 

- Review spectrometer data quickly on a CRT terminal . 

- Plot data on high resolution Varian printer/plotter in 
addition to line printer . 

- Plot, print or punch numerical agronomic and meteorological 
data stored on tape . 

- Fit polynomial curves of one to ten degrees to data and 
print regression coefficients and correlation coefficients. 

- Apply arithmetic transformations to data. 

- Cluster spectrometer data with one to 30 user specified 
bands and print tables with summaries of results. 

For the graphical software routines, EXOSYSDV makes use of the Graphics 
Compatibility System (GCS) acquired from the U.S. Army Corps of Engineers 
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at the Waterways Experiment Station, Vicksburg, Mississippi. There 
VN^as no cost for acquiring the GCS software package except for postage. 

The package of GCS software routines is very extensive for graphical 
purposes; the use of the GCS software routines to do the graphical 
plotting instead of writing completely new software reduced the imple- 
mentation time for the above mentioned tools from years to months. 

Examples of some of the graphical output are illustrated in Figures 
B-2 through B-5. Other graphical options available to the researcher 
using EXOSYSDV include polar coordinate system graphs, natural logarithmic 
and common logarithmic (base ten) formatted graphs, spline curve 
fitting, and other scaling and axes variations. 

The clustering capability implemented in EXOYSDV is the same 
cluster algorithm used in LARSYS (unsupervised with Euclidean distance 
measure) except that the Initial cluster centers are placed along the 
principal eigenvector. Examples of some of the new statistical output 
from EXOSYSDV are illustrated in Figures B-6 through B-S. Also the spectral 
information and some of the agronomic measurements stored on the 
spectrometer tapes can be punched and used in other statistical analysis 
packages such as SPSS or in user written Fortran programs. 

Documentation of the new capabilities in EXOSYSDV is the next step 
to be developed to made the system more useful to researchers both 
at Purdtie/LARS and to researchers at LARS' remote computer terminal 
sites. Initial development of documentation was completed this past 
year. In October, a gro,up of LARS researchers presented, in a one 
day session, an overview of EXOSYSDV to personnel at NASA/ JSC and 
demonstrated to several individuals the use of the system on the 
Purdue/LARS remote terminals at JSC. Development of new analysis 
tools and documentation x^ill continue so that researchers are able to 
make the fullest use of the data in the field research data library. 

In addition to the graphics and statistics software additions, the 
software for processing spectrometer data has been improved. The 
software for reformatting the Exotech 20C data xi'as revised to make it 
more efficient, saving both computer and personnel time. The original 



reformatting system was designed to handle a few hundred measurements 
per year, compared to 2,000 to 3,000 observations currently collected 
each year. At the same time additional items of agronomic and 
measurement data were added to the header record. Capabilities for 
updating the spectrometer data tapes have been added and documentation 
describing the reformatting systems written. The documentation will 
be very helpful in training new personnel how to use the systems. 


The design of the software system to 
radiometer and associated measurement and 
Programming will be completed next year, 
forerunner to the reformatting system for 
data (see Section C.2 ). 


process Exotech 100 
agronomic data was completed. 
This system will be a 
the new multiband radiometer 
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WAVELENGTH (mN) 


Figure B-2. Example graph from EXOSYSDV of spectral response as a 


function of wavelength. 



LEAF AREA INDEX 


Figure B^3. Example plot from EXOSYSDV of average spectral response for 
a user selected wavelength band versus an agronomic 
measurement using EXOSYSDV. 


Figure 


Figure 


1977 SPRING WHERT EXPERIMENT 



B-4. Example of one agronomic variable plotted against a 
second agronomic variable using EXOSYSDV. 


R2= 0.6243 Y= 23.916+ 4.541X1 



B-5. Example of polynomial curve fitting using EXOSYDV software. 
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Figure B-8L Example of the statistics printed for each cluster class using the new statistics 
capabilities in EXOSYSDV. 
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C. Multicrop Field Research 


5^7 

c-1 


Although the LACIE Field Measurements project (Bauer ^ al. , 1978)^ 
acquired a large quantity of spectral and agronomic data and analyses 
of it are providing quantitative information about the spectral 
characteristics of crops, the data are limited tp crops grovm in 
the Great Plains region of the United States. As we look ahead to 
the development of a global food and fiber information system utilizing 
remote sensing techniques, it is important to begin to conduct the field 
research required to more fully understand the spectral characteristics 
of crops other than wheat such as corn, soybeans, rice, cotton, and 
rangeland . 


As part of the LACIE Transition and Multicrop prbgrkm initiated 
in 1978, field research for com and soybeans was initiated. The 
overall objectives of the multicrop supporting field research are: 

— To increase the understanding of the basic radiation patterns 
of agricultural crops and their soil backgrounds under normal 
and stressed conditions. 

- To assess the capability of current and planned satellite 
sensor systems to capture available useful spectral infor- 
mation. 

- To determine the potential incremental improvement in 
performance of improved future systems. 


* . 

This section describing the results of work conducted as part of Task 2B, 
Application and Evaluation of Landsat Training, Classification, and Area 
Estimation Procedures for Crop Inventory, was prepared by C.S.T. Daughtry 
and B.F. Robinson. The Multicrop Field Research project was led by Dr. 
Marvin Bauer. The field measurements of corn and soybeans were directed by 
Dr. Craig Daughtry. Barrett Robinson and Larry Biehl were responsible for 
the spectral measurements. J4ff Kollenkark, Greg Walburg , David Longer, 
Michael Stabenfeldt, and Vic Fletcher assisted in the data collection. The 
field research instrumentation task was directed by Barrett Robinson. 
Professbrs David DeWitt and Leroy Silva made major contributions to develop- 
ment of the instrument system specifications. 

^Bauer, M.E., M.C. McEwen, W.A. Malila, and J.C. Harlan. 1978. Design, 
Implementation, and Results of LACIE Field Research. Proceedings , LACIE 
Symposium. October 23-26, 1978. NASA, Johnson Space Center, Houston, Texas. 
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The multicrop field research at Purdue/LARS during this contract 
year consisted of two components. The first was the acquisition of 
data describing the spectral properties of corn and soybeans under normal 
and stressed conditions. These experiments, conducted at the Purdue 
Agronomy Farm, are described in Section 1. 

The second component of the project, described in Section 2, was develop- 
ment of specifications for a new approach to spectral measurements for 
field research. The primary sensors used for LACIE Field Measurements 
were spectrometers capable of producing high-resolution spectra. 

In the future a new approach to the collection of field measurements data 
will be needed since it will not be feasible simply to multiply the 
current approach by the increased number of crops and regions which 
should be included in future experiments. Multiband radiometer systems 
can economically provide the necessary spectral measurements. With 
these instruments it will be possible to acquire measurements at 
more sites than is possible with the currently available high-spectral 
resolution spectrometer systems. And, it is more observations of 
crops and soils under a wide variety of conditions (not detailed 
spectral measurements of a limited number of locations and crop 
conditions) that is needed to increase our understanding of the spectral 
characteristics of agricultural scenes. There will be a continuing 
need for the high-resolution spectrometer systems to be utilized in 
field research, but less complex systems are also required which can be 
used to make large numbers of measurements at many sites economically and 
accurately. 



1. MEASUREMENTS OF SPECTRAL CHARACTERISTICS OF CORN AND SOYBEANS 


The potential of remote sensing to provide Information on the 
condition and yields of crops has not been adequately explored and 
developed, especially for corn and soybeans. This task represents 
the Initial phase of a multiyear research effort to evaluate and unaer- 
stand information about corn and soybeans crops contained in various 
types of remotely sensed data. The results of these carefully 
controlled field experiments can provide the foundation on which larger 
experiments Involving Landsat data should be built. 

1.1 Objectives 

The overall objectives of this field research task which include 
data acquisition and analysis are; 

1. To determine the reflectance characteristics of corn and 
soybeans as a function of maturity stage and amount of 
vegetation present. 

2. To examine the effects of moisture and nutrient stresses 
on the reflectance and radiant temperatures of corn and 
soybeans. 

3. To determine the effect of important cultural/management 
practices on the reflectance of soybeans. 

4. To develop and assess various methods of incorporating 
spectral, meteorological, and ancillary data into models 
of crop condition assessment and yield prediction. 

Activities this year towards fulfillment of these objectives 
include collection of spectral, agronomic, and meterological data at 
the Purdue University Agronomy Farm for five experiments with two 
different spectrometer systems. 

1.2 Description of Experiments and Measurements 

Several specific experiments have been established at the Purdue 
Agronomy Farm to address the above objectives. The experiments , 
described in the implementation plan, were: (1) corn moisture stress. 



(2) nitrogen fertilization of com, (3 and 4) phosphorous and potassium 
fertilization of corn and soybeans, and (5) cultural practices of 
soybeans. The experiments provide both normal and stressed conditions 
of corn and soybeans. Spectral and agronomic data were collected at 
approximately weekly intervals throughout the growing season. 

Spectral measurements were made with the Exotech 20C high 
resolution spectrometer system and the Exoteeh 100 Landsat-band 
radiometer (Figure C-1). The Exotech 20C field system collected data 
over the 0.4 to 2.4 pm range. The Exotech 100 field system which is 
a prototype of a more advanced multiband radiometer systems, (see 
Section 2), collected data in the four Landsat MSS bands .5~.6, 

.6-.7, .7-.8, .8-1.1 pm. The spectral data were calibrated using 1.2 
meter square painted barium sulfate standards. 

Augmenting the reflective measurements were radiant temperature 
measurements collected by Bames PRT-5 instruments and oblique- ground 
level and overhead photographs of the plots. The meteorological data 
included air temperature, barometric pressure, relative humidity, wind 
speed, and wind direction. A record of the irradiance was collected 
on pyranometer strip charts. Additional environmental data were 
acquired hourly by a computerized agricultural weather station located 
on the Agronomy Farm. 

Detailed agronomic measurements of the crop canopies included 
leaf area index, fresh biomass, dry biomass of leaves, stems, and ears 
(pods for soybeans), heights, maturity stage, percent of soil cover, and 
soil moisture. Leaf samples were collected and dried for nutrient 
analysis. Grain yields were measured at harvest. 

The approach to each of the objectives is briefly summarized 
below. The first of the overall objectives of this task was accomplished 
in the course of conducting specific experiments for Objectives 2 
and 3. Spectra of normal (unstressed) corn and soybeans will be assembled 
and plotted to, represent their reflectance as a function of maturity 
stage, biomass, If^af area index, height, and percent soil cover. 



Figure C-1. Spectroradlometer (top) and multlband radiometer (bottom) 

systems used for measuring spectral characteristics of 
corn and soybeans. 


To accomplish the second objective four experiments were conducted; 
one on corn moisture stress and three on nutrient stresses. Spectral 
and agronomic data were acquired at major maturity stages. The 
moisture stress experiment (Figure C- 1) was established on the Managed 
Soil Moisture System (MSMS). The MSMS facility is a series of 
plastic-lined plots equipped with pipes and pumps so that water and/or 
nutrient supply can be controlled and either type of stress can be 
imposed on the crop (Reetz, et al. 1978). Spectral reflectance, 
radiant temperatures, and agronomic characterizations of the stressed 
and nons tressed canopies were acquired once per day at seven to 10 
day Intervals. On selected days during the season, reflectance, radiant 
temperature, and leaf water potential were measured four to eight times 
to monitor the diurnal changes in these parameters due to stress . 

Nutrient deficiencies of corn and soybeans were studied on a 
Raub silt loam soil on which crop responses to residual and applied 
nutrients have been monitored for more than 20 years (Barber, 1970). 

Corn growth, development, grain yields and spectral reflectance were 
measured at four levels of applied nitrogen fertilizer (Figure C-2). 
Spectral reflectance of corn and soybean canopies growing at 
different levels of soil phosphorous and potassium also were examined 
at their major maturity stages from seedling to harvest (Figure C-3) . 

The third experiment examined the effects of different plant 
types, planting densities (within row spacings), and row spacing 
on the reflectance of soybeans (Figure C-4). Reflectance data acquired 
with the Exotech 100 system along with measurements of canopy variables 
(biomass, height, ground cover, etc.) will be used not only for character! 
zing the reflectance of soybeans, but also for examining the physical 
basis for increased soybean production at narrow row spacings. 

The fourth objective involved the assimilation of the data collected 
in these experiments into information which can be used in crop models 
to evaluate crop growth and to predict yields during the growing season. 
This objective represents one of the long range goals of this task and 
during the contract year consisted primarily of the preliminary conceptual 
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Figure C-3, Design and treatment descriptions of corn nitrogen fertilization 
experiment. Spectral measurements were made of pilots marked 
x^ith an asterik . 
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Figure C-4. Design and treatment descriptions of phosphorous and 

potassium fertilisation experiments . Spectral measurements 
were nuide on plots marked with an asterik 








Figure C-4. continued 
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Treatments are the entire amounts of fertilizer applied during four 
year rotation (com-soybeans-wheat-corn) . 
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Figure C-5, Design and treatment descriptions of soybean cultural practices 
experiment . 
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development of crop yield models which will incorporate spectral, 
meteorological, and ancillary data. 

Wlien the data acquired in this task are processed, various 
types of plots of the reflectance, radiant temperature, agronomic, 
and meteorological data will be made to verify data quality and to 
assess qualitatively the information contained in the data. Spectral 
data will be represented in several ways for analysis: band means 

for Landsat and the proposed thematic mapper bands; and transfor- 
mations of the reflectance values such as IR/red ratio and greenness- 
brightness. Canopy radiant temperature data will be represeiited as 
absolute temperatures, as relative temperature differences (canopy 
temperature minus air temperature), and as accumulated relative differences 
in temperature during critical growth phases. (Idso, et al. , 1977; 

Jackson et al. , 1977). Regression and correlation analyses will be 
used to relate biological and physical parameters of the canopies, 
such as leaf area index, biomass, height, leaf water potential, 
and maturity stage to measures of spectral response and radiant 
temperatures. Analyses of variance and covariance will be used to 
determine the threshold of detection and the separability of various 
levels of stress from each other. 

1.3 Data Collection and Processing 

Spectral measurements, along with agronomic and meteorological 
data, have been acquired on every day that weather conditions permitted. 

A summary of the data collection (number of observations per experiment) 
is shown in Tables C-1 and C-2. Crop maturity stages from seedling to 
senescence for 1978 are represented in these data. 

Thu cool wet weather of the spring of 1978 delayed soil preparation 
and corn planting at the Agronomy Farm and necessitated some modifications 
in the experiments planned. The Russell soil on which the corn 
irrigation experiment was to be planted remained too wet for plowing 
until earlyJune. On June 9 (a full month after the optimum planting 
data for corn) the first corn of this experiment was planted. This 
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date was too late to include a series of later planting dates as 
indicated in the implementation plan. Instead an experiment on the 
effects of soil surface moisture was conducted. 

Reflectance and radiant temperature data were collected over eight 
plots of corn at two different maturity stages with contrasting amounts 
of soil cover. 

More frequent rains then normal during the summer also made 
establishment and maintenance of different levels of water stress 
nearly Impossible. Short-term differences in soil moisture were 
established on several occasions in late July and August and radiant 
temperature and water potential were measured. However, much valuable 
experience on growing corn in the sand of the MSMS plots was gained , 
which should contribute more efficient operation of this facility in 
1979. 

The data acquired during the summer and fall are currently being 
processed and prepared for analysis. Analysis of the data is planned 
to start in early 1979. 

1.4 References 

Barber, S.A. 1970. Effect of rate and placement of phosphorous 

and potassium on a rotation, 18 year period. Research Progress Report 

365, Purdue University Agric. Expt . Station, West Lafayette, Indiana. 

Idso, S.B. R.D. Jackson, and R.J. Reginato, 1977. Remote sensing of 
crop yields. Science 196:19-25. 

Jackson, R.D. R.J. Reginato, and S.B. Idso 1977. Wheat canopy temperature: 
A practical tool for evaluating water requirements. Water Resources 
Research 13:651-656. 

Reetz, H.F. Hodges , and R.F. Dale. 1978. Managed soil moisture system for 
studying plant water relations under field conditions. Crop Science 
18: (in press) . 



Table C-1. Suramary of 1978 data acquisitioa on corn and soybean experiments 

(number of observations) by the Exotech 100 data acquisition system. 


Measurement 

Dates 


Corn 


Soybean 

Moisture 

Stress 

Phosphorous 

-Potassium 

Fertilization 

Nitrogen 

Fertilization 

Phosphorous 

-Potassium 

Fertilization 

Management 

Practices 

June 18-24 

96 




54 

June 25- July 1 

■ - 

92 

48 

66 

- 

July 2-8 

84 

70 

48 

84 

332 

July 9-15 

48 

46 

26 

64 

162 

July 16-22 

- 

- 

- 

- 

171 

July 23-29 

- 

- 

- 

- 

162 

July 30-Aug. 5 

- 

- 

- 

22 

218 

August 6-12 

- 

46 

26 

42 

54 

August 13-19 

- 

46 

20 

- 

108 

August 20-26 

16 

46 

- 

40 

162 

August 27-Sept. 2 

- 

- 

- 

- 

162 

September 3-9 

- 

46 

26 

42 

- 

September 10-17 

- 

- 

- 

- 

- 

September 17-23 

- 

- 


- 

108 

October 15-21 

' 


— 

— 

136 

Seasonal Totals 

244 

392 

246 

360 

1829 


Grand Total 3071 


Table C-2, Summary of data 1978 acquisition on com and soybean experiments 
by the Exotech 20C data acquisition system. 




Corn 



Soybean 


Dates 

Moisture 

Stress 

Phosphorous 

-Potassium 

Fertilization 

Nitrogen 

Fertilization 

Canopy 

-Soil 

Phosphorous 

-Potassium 

Fertilization 

Management 

Practices 

June 11-17 

28 

14 

- 

- 

— 

— 

June 18-24 

- 

- 

- 

- 

- 

- 

June 2 5- July 1 

6 

- 

20 

- 

10 

- 

July 2-8 

6 

48 

33 

- 

56 


July 9-15 

97 

16 

15 

- 

6 

- 

July 16-22 

8 

16 

- 

16 

29 

- 

July 23-29 

- 

18 

17 

- 

- 

- 

July 30-August 5 

8 

29 

17 

- 

25 

- 

August 6-12 

41 

- 

- 

- 

- 

- 

August 13-19 

— 

- 

16 

- 

5 

- 

August 20-26 

8 

32 

17 

16 

28 

- 

August 2 7- Sept. 2 

4 

- 

17 

- 

- 

- 

September 3-9 

. - 

32 

- 

- 

28 

- 

September 10-16 

- • 

- 

17 

- 

10 

- 

September 16-23 


— 

17 

— 

15 

33 

Seasonal Totals 

206 

205 

186 

32 

212 

33 


Grand Total 874 
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IP 


SPECIFICATION OF A STANDARDIZED MULTISPECTRAL DATA ACQUISITION 
SYSTEM FOR FIELD RESEARCH. 


2. 1 Introduction 

To develop the full potential of inultispectral data acquired from 
satellites, increased knowledge and understanding of the spectral 
characteristics of specific earth features is required. Knowledge of 
the relationships between the spectral characteristics and important 
parameters of earth surface features can best be obtained by carefully 
controlled studies over areas, fields, or plots where complete data 
describing the condition of targets is attainable and where frequent, 
timely spectral measurements can be obtained. The currently available 
instrumentation systems are either inadequate or too costly to obtain 
these data. Additionally, there is a critical need for standardized 
acquisition and calibration procedures to ensure the validity and 
comparability of data. 

The overall, long-term objective of this project is to develop a 
multispectral data acquisition system which will improve and advance 
the capability for remote sensing field research. The specific objectives 
are to: 

1. Specify, develop and test the prototype of a 
radiometric instrument system. 

2. Develop and document calibration, measurement and 
operation procedures. 

3. Develop software for data handling capability. 

The radiometric instrument will be a multiband radiometer with up 
to 8 bands between 0.4 and 2. 4 pm. The. data acquisition system will 
record data from the multiband radiometer, a radiation thermometer, and 
ancillary sources. The radiometer and data handling systems will be 
adaptable to helicopter, truck, or tripod platforms. 

The general characteristics of the system are that it will be: 

(1) comparatively inexpensive to acquire, maintain, and operate; 



(2) simple to operate and calibrate; (3) complete with data handling 
hardware and software; and (4) well-documented for use by researchers 

The specific results of the project will be: (1) a multiband 

radiometer; (2) a data acquisition and handling system; (3) a machine 
independent software package; and (4) a comprehensive system manual 
documenting the design and use of the above. The most significant 
result will be the establishment of the capability for researchers to 
acquire the data to effectively investigate relationships between the 
physical-biological and the multispectral characteristics of crops, 
soils, forests, water, and geological features. 

The instrument system will be a prototype of an economical, 
standardized system which can be utilised by many researchers to 
obtain large numbers of accurate, calibrated spectral measurements. 

As such it is a key element in Improving and advancing the capability 
for field research in remote sensing. 

2.2 General Characteristics of the Instrument System 

The complete Instrument system will provide the researcher with 
the means to acquire multispectral data with near optimum coverage 
of the spectrum from a variety of platforms ranging from hand-held 
operation to operation from a helicopter or small aircraft. The 
data playback feature will enable the researcher to check completely 
the performance of the system and measurement procedures in the field 
The data playback feature will also serve as a means to enter the 
gathered data directly into a digital computer using the acquisition 
module as the interface. 

Specifications for the instrument system will be established 
with adherence to four principles to enhance the eventual commercial 
availability of the system and adoption of the system by the remote 
sensing community. 



(1) The specifications will emphasize a modular design. 

The design will allow for system expansion or alteration 
through addition or substitution of modular building blocks. 
Design requirements will Include that the Instruments be 
functionally flexible and easy to understand, trouble-shoot, 
and repair. 

(2) The final system will be assembled from commercially available 
(from the chosen vendors) modules and will require 
essentially no electronic or mechanical construction. 

(3) Only hardware employing accepted signal processing 
techniques will be specified in the system. 

(4) The final system must be easily operated by pei’sons 
with diverse backgrounds and limited electronics 
training and easily repaired by persons with minimal 
electronics training. 


The modular design approach will make possible the purchase of the 
particular system best matched to individual research goals and resources. 
The lowest cost (and least flexible) system would require only the 
radiometer module plus a multiple position switch and digital voltmeter. 
With this simple system, radiometric data would be acquired by sampling 
each output channel of the radiometer and manually recording the voltage. 

A second, more complete system would require the radiometer module, 
plus a data handling module. The data handling module would convert 
the analog signals from the radiometer to digital numbers and print them 
on paper tape for later analysis. 

A third system would entail expansion of the data system to 
include a digital recording and playback capability. In this 
expanded system the paper tape data optionally printed by the data 
handling module would serve to verify (in the field) the accuracy of 
the data stored by the data storage medium. 

The general specifications for each of the modules follow: 


Multiband Radiometer Module 


The purpose of the multiband radiometer module is to simultaneously 
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produce analog voltages which are proportional to the scene radiance 
in each of the spectral bands. The general specifications of the radiometer 
module are : 


- Rugged, lightweight, and portable 

- Battery operated 

- Relatively low cost 

- Eight spectral bands from 0.4 to 12.5 micrometers 
(Tentative band selection is 1.45-0.52, 0.52-0.58, 
0.63-0.69, 0.76-0.90, 1.10-1.30, 1.55-1.75, 2.08-2.35, 
10.5-12.5. This selection includes the seven Thematic 
Mapper bands and an additional band in the near infrared.) 

- Optical filters to define each spectral band 

- Optical units to be exchangeable to obtain the spectral 
pass bands desired by the researcher 

- Selectable field of view (by lens substitution) 

- All radiometric data simultaneously available in parallel 
format. 


Data Handling Module 

The data handling module will convert the signals from the radiometer 
to digital format and store them with digital ancillary data (i.e., 
date, time, observation number) in a removeable solid state module 
having a data retention battery. The data handling module will print 

data, compute reflectance or radiance, and also serve as the interface 
to enter data to digital computers and peripherals. The general 
specifications of the data handling module are: 

- Rugged, lightweight, and portable 

- Battery operated 

- Relatively low cost 

- Parallel multichannel analog input port compatible with analog 
outputs of radiometer 

1 

- A simultaneous sample and hold circuit for each detector 
channel (or equivalent) 

- Remote activation of the acquisition sequence enabling 
synchronous operation with boresighted cameras operated by 
intervalometers 

- Provision to enter date, time and observation number 
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- Digital hard copy printer to tabulate data values of each 
radiometer channel, plus ancillary data 

- Digital data output port for RS 232-C and "12-blt parallel" 
with request to send/clear to send lines 

- Auxllllary analog input ports to accomodate signals from 
other sensors such as the Barnes PRT-5 thermal sensor 

- Provision for reflectance factor or radiance computations. 

Documentation and User Training 

In order for researchers to acquire accurate and comparable 
measurements it will be important to fully document and describe use 
of the instrument system. An educational package including audio- 
tutorial materials is being planned. It would provide researchers with 
the theoretical and practical knowledge to gather, process, and analyze 
high quality spectral data. Major topics to be treated are: 

(1) Measurement and calibration fundamentals and procedures, 
including important parameters of the measurement 
situation and procedures for data evaluation and comparison. 

(2) Specifications for the radiometer and instrument system, 
system tests, and typical results. 

(3) Methods of positioning the instrument, including suggested 
hardware and mounting techniques for operation on a 
helicopter, mobile tower, and tripod platforms. 

(4) Data handling and data processing procedures including 
a machine independent computer program to organize and 
present the data. 

(5) Case studies of typical experiments including experimental 
design, measurement and data handling procedures, and 
analysis techniques and results. 

2.3 Objectives 

The objective for this contract year was to coordinate development 
of the specifications and the preparation of a draft request for 
quotation (RFQ) for a portable. Inexpensive field Instrument system 
with multiple, selectable, spectral bands to measure, record, and 
organize large amounts of high quality spectral data. 


2.4 Approach 


This effort I'as the initial step required to acquire 
instrumentation needed for meaningful field research; i.e. , develop- 
ment of an inexpensive standardized, versatile field instrument to 
provide calibrated intercomparable spectral data when used by researchers 
at different sites. 

The general characteristics of the radiometer and data acquisition 
system were presented to and discussed with active remote sensing 
researchers at other institutions with regard to their experimental 
needs and data handling preferences and capabilities. Potential 
vendors were contacted for their reaction to the specifications. Vendors 
were encouraged to suggest alternate means to obtain a practical, 
effective, commercially available hardware system. Final specifications 
will be prepared in cooperation with personnel from NASA/ JSC. 

A fonnal request for quotation will be prepared for submission 
to qualified vendors. The goal will be to ensure that each subsystem 
will be available for purchase on a continuing basis from a reliable 
supplier. The document will be prepared in cooperation with personnel 
from NASA/ JSC. 

Initial Discussions 

Initial discussions held with Earth Observation Division (EOD) and 
Earth Resources Program Office administrators at NASA/JSC in the 
Spring of 1978 indicated that the proposed system could play a major role 
in supporting future research in remote sensing. As a result of these 
discussions, the EOD directed LARS to undertake, in the summer and fall 
of 1978, the task of developing the specifications of the instrument to 
the point where a formal RFQ can be prepared. 

Contacts with Researchers and Potential Users 

On May 24, 1978, a meeting of NASA and USDA researchers was hosted 
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by the USDA/ Science and Education Administration facility at Weslaco, 
Texas. Discussions of the preliminary specifications were led by 
Purdue/LARS researchers and the need for the system was confirmed. 

During the contract year, the specifications were also discussed with 
other scientists from England, France, Germany, and the U.S. Their 
comments and suggestions were considered in the definition of the speci- 
fications. In each case, the need for the system was expressed by the 
researcher. 

Following the contacts with researchers and vendors, a more 
detailed set of preliminary specifications and system configuration 
alternatives was prepared and presented to the participants of the 
Weslaco meeting and to researchers from ERIM and the University of 
Nebraska, On October 12, a meeting was held at Purdue/LARS to discuss 
these specifications. 

Contacts with Potential Vendors 

In June 1978 potential vendors were provided copies of an 
"Informal Request for Consideration" which listed potential specifications 
for the Data Logger and the Multiband Radiometer Modules. During the 
summer, Professors Silva and DeWitt made follow-up calls to interested 
vendors. Professor Silva visited eight vendors to discuss the data 
logger. Professor DeWitt visited two vendors v^ith regard to the 
multiband radiometer. A third Vendor interested in, developing the 
multiband radiometer module has been in telephone contact. 

Preparation of Specifications 

Preparation of tentative detailed specifications was completed in 
August and the specifications for the Radiometer and Data Module were 
fot'warded to interested manufacturers and potential users for their 
reactions and suggestions. The first draft of the specifications for 
the RFQ was completed in September 1978. 


Significant developments in the specifications for the multiband 
radiometer are the probable inclusion of a thermal infrared channel, 

10.5 to 12.5 ym, and the use of temperature compensated detectors. 
Significant developments in the specifications for the data logger 
are the capacity for direct reflectance computation and the use of 
removeable semi-conductor data storage (instead of magnetic tape). 

As veil as providing more easily interpreted feedback in the field, 
these features will enable many researchers to efficiently process 
all their data to bidirectional in-band reflectance factor (or thermal 
band radiance) using only the data acquisition module. 

Field Tests 

An experimental truck-mounted boom vas constructed and tested 
for measurement of Landsat band bidirecticmal reflectance factor 
and radiant temperature. An Exotech Model 100, a Barnes PRT-5 and 
a motor dirve camera were mounted on the boom (see Figure C-1) . The 
system was used for measurements of corn and soybean canopies from 
altitudes of 6 and A meters, respectively. The routine use of this 
system provided experience which will lead to improved design for an 
easy to construct, inexpensive, pick-up truck operated system which is 
sufficiently versatile to serve a large number of researchers. 

Helicopter operation of the Erotech Model 100 boresighted with 
a camera was accomplished with the funds of other agencies and projects. 
Calibration was accomplished on-site with a canvas reflectance panel. 

The instrument and data acquisition system performed correctly and it 
was demonstrated that an untrained crew could quickly learn to perform 
the calibration and data acquisition operations. 



2.5 Plans for Next Contract Year 


During the next year the following tasks and accomplishments 
leading to commercial availability and use of the instruments in 1980 
are planned; 

A review panel will be established to review the 
development, procurement and testing of the system 

A prototype system will be acquired and tested by 
Purdue/LARS 

- Based on tests of the prototype system, final specifications 
for production units will be prepared 

A system manual will be prepared 

An experimenters’ manual will be prepared 

Software for storage and use of radiometric and ancillary 
data will be developed. 


D. Determining the Climatic and Genetic Effects 
on the Relationships between Multispectral Reflectance 

A 

and Physical-Chemical Properties of Soils 

1. INTRODUCTION 

Although a large body of knowledge has been accumulated about the 
physical and chemical characteristics of soils as they are Influenced 
by the soil forming ^actors of climate, parent material, relief, biologi- 
cal activity, and time, there is only limited knot<7ledge of how these 
factors relate to the reflected radiation from surface soils. Earlier 
studies have shovm that information about the spectral properties of 
soils may be useful in their identification and characterization (1,2, 3, 4). 

If present satellite sensors and the improved sensor systems 
planned for future satellites are to be used most effectively in the pre- 
paration of land use capability maps and soil productivity ratings as 
these relate to crop production, it is crucial to define quantitatively 
the soil variables related to productivity which can be measured by or 
correlated with multispectral radiation from the surface soil. This 
research seeks to contribute significantly to the understanding of the 
multispectral reflectance of soils as it relates to climate, physical 
and chemical properties of soils , engineering aspects of soil use and 
potential agricultural productivity. 


This section describing the results of Task 1C, Determining the Climatic 
and Genetic Effects of the Relationships between Multispectral Reflectance 
and Thysieal-Ghetnieal Properties of Soils, was prepared by the task manager 
Eric Stoner. Prof. Marion Baumgardner and Dr. Richard Weismiller are co- 
leaders of this task. Barrett Robinson and Larry Biehl were responsible 
for the spectral measurements. Lou Nash, Steve Jordan, Charles Baker, 

Lyn Kirschner, and staff of the Purdue Soil Characterisation Lab made the 
chemical and physical measurements. Prof. Virgil Anderson assisted with 
the experiment design. 


1.1 Quantification of Soil Properties 

Modern soil classification systems emphasize the importance of infor- 
mation about the quantitative compositions of soils. In order to differ- 
entiate among soil groups, it is necessary to rely on laboratory measure- 
ments of selected soil properties. Physical, chemical, and engineering 
determinations of most soil properties follow well established procedures 
of laboratory analyses. Certain of these soil properties are selected as 

diagnostic criteria in the soil classification process, based on their 
Importance in understanding the genesis of the soil. By a procedure of 
empirical correlation, critical limits between sets of soils are estab- 
lished, designed to reflect the influence of the soil forming factors of 
climate, parent material, relief, biological activity, and time. 

Quantitative measurements of soil spectral properties have become 
available as a diagnostic tool for the soil scientist with the advent of 
such instruments as the Exotech Model 20C spectroradiometer . However, the 
climatic and genetic effects on the relationships between measured spec- 
tral properties and specific chemical, physical, and biological properties 
of the soil are not well understood, l^ereas soil color is used as diag- 
nostic criterion in the U.S. Soil Taxonomy (5), the determination of soil 
color by comparison with a color chart continues to be a rather nonquan- 
titative and subjective procedure. Spectral characterization of soil 
"color" by means of quantitative spectroradiometric measurements may add 
to the precision with which soils can be differentiated. With this in-^ 
creased precision of soil spectral characterization, the relationships 
with the more important diagnostic soil characteristics or qualities that 
are not so easily and accurately observed may be better understood. 

2. STUDY OBJECTIVES 

The general objective is to define quantitatively the relationships 
between soil reflectance and physicochemical properties of soils of sig- 
nificance to agriculture and engineering. Selection of soil samples with 
a wide range of important soil characteristics by statistical stratifica- 
tion of continental United States climatic zones permits the evaluation 
of climatic and genetic effects on the relationships between multispec- 
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tral reflectance and these soil properties. A further objective is to 
define the relationships sufficiently to design further research to 
quantify the contributions which different soil components make to the 
multispectral characteristics of specific soils. The ultimate objective 
of this research approach is to provide a body of knowledge and interpretive 
skills which will render remote multispectral sensing a valuable tool 
for mapping soils, determining land use capabilities and soil productivity 
ratings, identifying crops and predicting crop yields. 

3. EXPERIMENTAL APPROACH 

3.1 Stratification and Sampling 

Approximately 250 soils, representing a statistical sampling of the 
more than 10,000 soil series in the United States were selected for this 
investigation. Selections were made from a list of the more than 1300 
Benchmark soil series representing those soils with a large geographic 
extent and whose broad range of characteristics renders these soils so 
widely applicable for study. 

Stratification of soil sampling was based on series type location 
within climatic zones. Climatic strata included the frigid, mesic, thermic, 
and hyperthermic soil temperature regimes as defined by the U.S. Soil Taxon- 
omy (5,6,7) as well as the perhumid, humid, subhumid, semiarid, and arid 
moisture regions as identified by Thornthwaite ' s 1948 Moisture Index (8) , 

A random selection procedure was used within each stratified climatic 
zone to select a number of soils series approximately in proportion to 
the geographic extent of that region. Resulting sample distribution by 
climatic region for the soils actually received is shovni in Table D-1. 

Considerations*were also made to include soils which represent the 
major parent material categories and the ten soil orders of the U,S. 

Soil Taxonomy (5). Table D-2 presents the distribution of tie Benchmark 
soil series on-hand according to soil parent material. As can be seen in 
Table D-3, the distribution of Benchmark soils used for this study is very 
similar to the areal extent of the nine soil orders found in the continental 
United States (Oxisols being absent in the contiguous states). 
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The type of clay mineralogy present in a soil is also known to influ- 
ence its reflectance characteristics (3,4). The mineralogy classes, one 
of the family differentiae for mineral soils in the U.S. Soil Taxonomy (5), 
are listed for the Benchmark soils on-hand in Table D-4. 


Table D-1. Distribution of soils by climatic region. 


Number of Benchmark 

Climatic Region Soil Series 


1. 

Perhumic Mesic 

6 

2. 

Humid Frigid 

18 

3. 

Humid Mesic 

37 

4. 

Humid Thermic 

30 

5. 

Humid Hyperthermic 

6 

6 . 

Subhumid Frigid 

21 

7. 

Subhumid Mesic 

23 

8. 

Subhumid Thermic 

18 

9. 

Subhumid Hyperthermic 

2 

10, 

Semiarid Frigid 

9 

11. 

Semiarid Mesic 

24 

12. 

Semiarid Thermic 

10 

13. 

Semiarid Hyperthermic 

5 

14. 

Arid Frigid 

2 

15. 

Arid Mesic 

16 

16. 

Arid Thermic 

12 

17. 

Arid Hyperthermic 

1 


240 total 



Table D-2. Distribution of soils by 


Parent Material 


Rocks weathered in place 
Igneous 
Sedimentary 
Metamorphic 
Soft rock residuum 

Transported materials 

Alluvium (general) 

Calcareous 
Non-calcareous 
Colluvium 
Lacustrine 
Marine sediments 
Loess 

Other eolian sediments 
Glacial drift 
Till 

Calcareous till 
Glaciofluvial deposits 
Glacial outwash 
Glaclolacustrine materials 

Organic materials 
Loamy sediments 
Silty sediments 
Calcareous silt loam 
Marsh deposits 
Pedisediments 





Number of 

Benchmark Soil Series 


10 

27 

1 

5 


38 

11 

15 

3 

4 

16 
29 

8 

29 

5 
4 
4 
2 

2 

17 

4 

3 

1 

2 

240 total 
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Table D-3. Distribution of soils by Soil Order. 



Benchmark 

Soils 

United States Extent 


Number 

Percent 

Percent 

Mo Hi sol 

73 

30.4 

24.6 

Alfisol 

40 

16.7 

13.4 

Entisol 

39 

16.2 

7.9 

Aridisol 

25 

10.4 

11.5 

Ultisol 

22 

9.2 

12.9 

Inceptisol 

18 

7.5 

18.2 

Spodosol 

15 

6.2 

5.1 

Vertisol 

4 

1.7 

1.0 

Histosol 

4 

total 240 

1.7 

0.5 


Table D-4. Distribution of soils by mineralogy classification. 


Class 

Number of 

Benchmark Soil Series 

Mixed* 

176 

Montmorillonitic 

35 

Siliceous 

17 

Carbonatic 

5 

Kaolinitic 

3 

Gypsic 

2 

Micaceous 

2 


’"'Designates class with <40% of any one mineral other 
than quartz or feldspars. 
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3.2 Acquisition of Soil Samples 

The Soil Survey Investigations Division of the Soil Conservation 
Service (USDA) cooperated with LARS in the collection of field samples 
from 39 states. Duplicate field samples were collected for all Benchmark 
soil series requested; one sample from a site near the type location for 
the current official series, and one sample from a site located from one 
to 32 kilometers from the first site and in a different mapping delinea- 
tion. Soil Conservation Service field survey personnel were responsible 
for sample collection of Benchmark soils in their locality. Of the 
original list of approximately 250 Benchmark soils requested, the Soil 
Conservation Service has collected, properly identified, and forwarded 
240 Benchmark soils, or 480 duplicate soil samples to LARS. This excel- 
lent response of over 95 percent of the requested samples forms an out- 
standing collection of soil samples for detailed chemical, physical, and 
spectral analysis. All samples conform to the central concept of each 
individual soil series as each soil would be identified and mapped by an 
experienced soil surveyor in the field. 

3.3 Preparation of Soils for Analysis 

After receipt of the soil samples and initial data logging, samples 
were dried, crushed, and sieved to remove all particles larger than 2 ram 
diameter. Cardboard containers were used to store subsaraples of each 
soil sample for chemical, physical, spectral, and engineering determina- 
tions . 

3.4 Spectral Measurements 

The Exotech Model 2QC was used in an indoor configuration with the 
bidirectional reflectance factor ref lectoraeter (9, 10) in order to obtain 
spectral readings in the 0. 52-2. 32vmi wavelength range. The illumination 
source was a 1000 watt tungsten iodine coiled filament lamp which trans- 
fers a highly collimated beam by means of a paraboloidal mirror to the 
sample-viewing plane. Detector height above the sample was 2.44 m, and 
a 3/4° field of view required that the sample holder be approximately 
10 cm in diameter. 
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Because of the need to provide an equlpotential condition for spectro- 
radiometric analysis of the prepared soil samples, a procedure was chosen 
which creates a PF2 soil moisture tension on all the soil samples (11,12). 
Two asbestos tension tables were constructed and a 100 cm column of water 
was established to create a PF2 (approximately one-tenth bar) soil moisture 
tension for up to 56 soil samples at one time. Sample holders were de- 
signed and constructed of plastic rings 2 cm deep by 10 cm in diameter with 
60 mesh brass strainer cloth stretched taut and fastened in a countersunk 
groove in one end. Sample holders were painted with non-reflecting black 
paint to reduce unwanted reflection external to the target of interest. 

After saturation of the soil filled, level sample holders for about 
four hours, the samples were placed on the tension tables for 24 hours in 
order to reach equilibrium. The PF2 moisture tension was desirable mainly 
for the ease with which large numbers of samples could be prepared at equi- 
potential moisture characteristics. Shortly after placement of each sample 
holder on the sample table of the ref lactometer for spectral readings, a 
portion of the sample was transferred to a moisture tin, weighed, dried in 
a forced air oven at 105°C, weighed again, and moisture content reported 
as percentage of oven dry weight. The setup for indoor spectral measure- 
ments of soils is shown in Figure D-1. 

3.5 Measvirement of Physicochemical Properties 

Facilities of the Purdue University Agronomy Department were used to 
perform the following analyses: mechanical analysis (eight particle size 

separates), organic matter content, cation exchange capacity, and iron 
oxide content. Other associated measurements to be taken are shown in 
Figure D-2. Because of the late arrival of many of the soil samples, me- 
chanical analysis and iron oxide analysis are still in progress, while 
cation exchange capacity measurements have been completed. 

3. 6 Measurement of Engineering Properties 

The following engineering properties are being determined in the 
facilities of the Purdue University Department of Geosciences: liquid 

limit, plastic limit, plasticity index, activity, liquidity index, volu- 
metric shrinkage, and compression index. Coding of these measurements 
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b. Saturated sample being placed on asbestos t.insion table 


Figure D-1. Setup for laboratory spectral measuremen s of soils. 
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and resulting engineering soils classifications can be seen in Figure D-2. 

Of the 226 soils processed thus far, 173 or 77% have contained enough 
of the sub-. 425 mm fraction to run the battery of Atterberg and associated 
tests (13). The remainder of the soils are sands and exhibit virtually 
no cohesive properties at all. 

The measured data values for the engineering properties are represen- 
tative of the data values given on the established series description sheets 
for the soil series and agree well, where comparison is possible. Major 
variations in properties did occur in a few cases between the duplicate 
samples of a given soil, but in general, measurements of engineering proper- 
ties between duplicate samples were in agreement. 

3.7 Establishment of Data Logging Procedure 

An Identification record containing complete soil taxonomic information 
along with site characteristics and laboratory analysis results has been pre 
pared and implemented for storage and rapid retrieval as part of the EXOSYS 
software package (14). The record consists of seven computer cards of infer 
mation for each observation (Figure D-2) . 

4 . RESULTS OF SOIL SPECTRAL MEASUREMENTS 

The bulk of the soil spectral measurements were made on nine days from 
August 2,9 to September 15, 1978 with 56 soil samples being analyzed each 
day. Additional measurements were taken on October 20, 1978 for late- 
arriving soil samples. During this time period the Exotech Model 20C spec- 
troradlometer was being used for outdoor spectral measurements, but it was 
possible to dismount the detector head and position it in the optics labora- 
tory in the evening after taking outdoor daytime measurements. Connector 
cables were passed through a wall port for connection of the detector head 
with the recording instruments in a van parked outside. Careful alignment 
of the reflectometer setup with the detector helped to insure repeatability 
of experimental conditions. 

Reformatting of the spectral measurement data and filing of the data 
on a sequential run tape was completed on November 10, 1978. Preliminary 
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Figure D-2. EXOSYS soil record sheets for header record data logging 
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^analysis of the data has shown all measurements to be of excellent quality. 
' Initial data analysis included primarily graphical display of soil spectral 
measurements of Interest * Statistical analysis will follow upon comple- 
tion of header record information for each soil sample. 

4.1 Evaluation of Experimental Procedures 

During each of the ten days in which soil spectral measurements 
were taken, two check samples were randomly assigned to each block of 56 
soil samples as a verification of repeatability of the tension table setup 
as well as the instrument setup. The check samples consisted of subsamples 
of a larger field sample of Fincastle silt loam soil (Aerie Ochraqualf) . 
Figure D-3 presents the spectral curves for these 20 check samples, divided 
into two time blocks for sake of illustration. Soil moisture content, or 
more correctly, soil moisture weight percent (MW%) as measured gravimetri- 
cally after taking spectral readings is shown. As can be easily seen, 
moisture weight percent varied little from an average 31.3 MW% for all 
check samples, while the spectral response curves are remarkably similar 
for all days and setups. 

4.2 Representative Soils from Different Soil Orders, Climatic Zones, 
and Having Different Mineralogy Classes 

The diversity of soil spectral response is evident from the soil 
curves presented in Figure D-4. As listed in Table D-5 , all ten soil 
orders of the U.S. Soil Taxonomy are represented here, including four 
Oxisols from Brazil which were included to contrast with what has conven- 
tionally been thought of as the "typical” soil reflectance curve. 

Each pair of curves for a given soil series represents the duplicate 
£ie].d samples which may vary in certain soil characteristics within an 
allowable range permitted for that soil series. For example, small dif- 
ferences in clay content or organic matter Content may in themselves 
affect soil reflectance while at the same time directly influencing the 
moisture weight percentage in soil samples equilibrated at an equlpotential 
moisture tension. It is generally recognized that increased moisture 
weight percentage results in decreased soil reflectance throughout the 
0.4 to 2.5 urn wavelength range (1,4,15). This soil moisture influence 
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Figure D-3. Soil spectral curves and moisture weight percentages 
(MW%) for 20 check samples from ten different setups 
of the tension table apparatus with resultant reflec- 
tance measurements . 
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Figure D-4 . (Continued) 
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Figure D-4. (Continued) 
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Figure D-4. (Continued) 



Table I>-5. Representative soils of ten Soil Orders with various mineralogy classes and from different 
climatic zones (refer to Figure D-2) . 


Soil Series Name (State) 

Taxonomic Classification 
Order Subgroup 

Climatic Zone 

Mineralogy Class 

Superstition (AZ) 

Aridisol 

Typic Calciorthid 

Arid Hyperthermic 

Mixed 

Svea (ND) 

Mollisol 

Pachic Udic Haploboroll 

Subhumid Frigid 

Jlixed 

Marias (MT) 

Entisol 

Ustertic Torriorthent 

Semiarid Frigid 

Montmorillonitic 

Reno (NV) 

Aridisol 

Abruptic Xerollic Durargid 

Semiarid Mesic 

Montraorillonitic 

Newtonia (OK) 

Mollisol 

Typic Paleudoll 

Humid Thermic 

Mixed 

Victoria (TX) 

Vertisol 

. Udic Peliustert 

Subhumid Hyperthermic 

Montmorillonitic 

Pompa'no (FL) 

Entisol 

Typic Psaramaquent 

Humid Hyperthermic 

Siliceous 

Terra Ceia (FL) 

Histosol 

Typic Medisaprist 

Humid Hyperthermic 

(organic soil) 

Russell (IN) 

Alfisol 

Typic Hapludalf 

Humid Mesic 

Mixed 

Pacolet (SC) 

Ultisol 

Typic Hapludult 

Humid Thermic 

Kaolinitic 

Norwich (NY) 

Inceptisol 

Typic Fragiaquept 

Humid Mesic 

Mixed 

Pence (WI) 

Spodosol 

Typic Haplorthod 

Humid Frigid 

Mixed 

Cascavel (PR, Brazil) 

Oxisol 

Haplic Acrorthox 

Humid Thermic 

Oxidlc 

Pa to Branco (PR, Brazil) 

Oxisol 

Haplic Acrorthox 

Humid Thermic 

Kaolinitic 

Guarapuava (PR, Brazil) 

Oxisol 

Typic Acrohumox 

Humid Thermic 

Oxidic 

Londrlna (PR, Brazil) 

Oxisol 

Typic Haplorthox 

Humid Hyperthermic 

Kaolinitic 
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can be seen in many of the duplicate samples In Figure D-4, although moisture 
alone is not responsible for this overall ‘'darkening" effect. 

Aside from the obvious soil water absorption bands at 1.44 and 1.94um, 
a small, yet distinct absorption band can be seen at 2.2 um in the Pacolet 
(SC), Marias (MT) , and other soils. This band has been attributed to the 
hydroxyl group and is especially prominent in kaolinlte and montmorillonite 
clays. Another broad, but well defined absorption band is that at 0.9 um, 
attributed to iron oxides in the soil. This band is obvious in all four 
Oxlsols from Parana State, Brazil, some of which have iron oxide contents 
as high as 25%. 

4.3 Comparison of Laboratory and Field Spectral Data 

Results illustrated in Figure D-5 indicate that a comparison can be 
made between field measured Soil plots under solar illumination and labora- 
tory measured soil samples under artificial illumination when all spectral 
measurements are compared to a BaSO^ standard. The general trends of the 
curves for each soil are similar no matter whether the measurements were 
taken in the lab, or in the field under four different surface conditions: 
bare dry soil, bare moist soil, residue-covered dry soil, or residue-covered 
moist soil (cut-offs in field measured curves represent attenuation of solar 
illumination source) . This concept is an important one — that a piven soil 
has a wide range of spectral reflectance characteristics defined by specified 
conditions of irradiation and viewing and varying according to naturally 
occuring states of surface roughness, residue amount, and green vegetative 
cover. 

4.4 Comparison of Soils with Differing Engineering Properties 

Illustrations of the spectral curves of four soils of interest are 

given in Figure D-6. The pertinent behavioral data are given in Table D-6 . 
The Abbott IS a fine, montmorillonitic clay from the arid mesic climatic 
zone of Millard County, Utah. it experienced a volumetric shrinkage 
rate of 101 making it a soil which would be of some concern in foundation 
design due to the potential for a large amount of swelling upon wetting. 

This would be of particular concern in using slab foundations in this arid 


area. 
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FINCflSTLE SIL (flERIC OCHRfiQURLF ) 
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Figure D-5. Compatlson of two soils measured in the field and in the labora- 
tory. Percentage figures are moisture weight percent; RES = 
corn residue cover of 35%; BARE = residue free field plots; 

LAB = laboratory data. 
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Figure D-^-6 i Comparison of four soils with differing engineering properties. 



Property 

Abbott 

Dia 


Brackett 

Liquid limit 

61 

52 

39 

31 

31 

39 

Plastic limit 

25 

24 

22 

23 

23 

31 

Plasticity index 

36 

28 

17 

8 

8 

8 

Activity 

N/A 

N/A 

0.95 

0.50 

0.38 

0.28 

Shrinkage limit 

17 

15 

20 

19 

16 

23 

Shrinkage ratio 

1.8 

1.8 

1.7 

1.7 

1.8 

1.6 

Volumetric shrinkage 

101.0 

85.3 

42.0 

29.2 

41.4 

38.4 

Linear shrinkage 

20.1 

18.5 

11.0 

8.1 

10.8 

10.2 

Compression index 

.439 

.318 

.261 

.189 

.189 

.261 

ASTM Particle Size Distribution 
medium sand (2-. 425 mm) 

00.8 

00.9 

8.5 

5.4 

9.1 

8.2 

fine sand (.425-.075 mm) 

09.5 

10.1 

42.1 

54.3 

32.3 

19.1 

fines (<.075 mm) 

89.7 

89.0 

49.4 

40.3 

58.6 

72.7 

Specific gravity 

2.61 

2.43 

2.47 

2.48 

2.50 

2.44 

MSHO classification 

A-7 

A-7 

A-6 

A-6 

A-6 

A-6 

Unified classification 

CH 

CH 

CL 

ML 

ML 

ML 

Cation exchange capacity 

49.8 

44.4 

23.1 

26.7 

23.7 

26.7 

Clay content 

N/A 

N/A 

17,9 

16.1 

21 

L 

28.1 
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Dia, a fine loam from the arid raesic climatic zone of Churchill County, 
Nevada, and Brackett, a loam from the suhhumid thermic climatic zone of Bell 
County, Texas typify the majority of the soils studied with respect to 
their engineering properties. They exhibit moderate liquid limits and vol- 
umetric shrinkage readings, and are well-graded. Variation between 
duplicate samples are low. The spectral curves agree with what was ex- 
pected given the moisture content at which the readings were taken and the 
clay content of the sample. 

The Fortwingate soil is a friable loam from the semiarid mesic climatic 
zone of McKinley County, New Mexico. It exhibited the lowest volumetric 
shrinkage reading to date. It is better graded than the Abbott and has a 
lower cation exchange capacity. In comparing the Abbott soil with the Fort- 
wingate soil it is evident that the low shrink-swell soil, the Fortwingate, 
does not have as large a water absorption band at 1.4 urn as the high shrink- 
swell soil, the Abbott. This may indicate a spectral region of the soil 
response curve useful for the prediction of high shrink-swell conditions 
which would be of concern to the engineer . 

5. STATISTICAL ANALYSES 

Upon completion of the measurements of spectral, chemical, and 
physical properties, all data will be entered into the soils data record 
in preparation for the statistical analyses. Several approaches are 
anticipated: (1) analysis of the variability of individual soil properties 

relative to the duplicate samples of a given soil series, (2) analysis of 
the degree to which soil series can be uniquely distinguished for individual 
properties, and (3) analysis of the possiblity of separating Benchmark 
soils on the basis of their reflectance characteristics. 

Measured soil properties will be analyzed for normality and homogeneity 
of variances and appropiate transformations will be performed when needed. 
Analysis of variance will be used to test for significance of spectral 
differences between duplicate soil samples and among soil series. If war- 
ranted by the analysis of variance tests, multivariate analysis will be per- 
formed for relating the response in selected combinations of wavelength 
bands to Chemical, physical, and engineering properties for soil samples 
grouped according to climatic zone and parent material categories. 
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